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We systematically examine various proposals which aim at increasing the accu- 



o 
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■ racy in the determination of the renormalization of two-fermion lattice operators. We 

concentrate on three finite quantities which are particularly suitable for our study: 
the renormalization constants of the vector and axial currents and the ratio of the 
^ • renormalization constants of the scalar and pseudoscalar densities. We calculate these 

D , quantities in boosted perturbation theory, with several running boosted couplings, at 

' the "optimal" scale q* . We find that the results of boosted perturbation theory are 

usually (but not always) in better agreement with non-perturbative determinations of 
^ , the renormalization constants than those obtained with standard perturbation theory. 

H ' The finite renormalization constants of two-fermion lattice operators are also obtained 

non-perturbatively, using Ward Identities, both with the Wilson and the tree-level 
Clover improved actions, at fixed cutoff {(3 = 6.4 and 6.0 respectively). In order to 
amplify finite cutoff effects, the quark masses (in lattice units) are varied in a large in- 
terval ^ am ^ 1. We find that discretization effects are always large with the Wilson 
action, despite our relatively small value of the lattice spacing (a~^ ~ 3.7 GeV). With 
the Clover action discretization errors are significantly reduced at small quark mass, 
even though our lattice spacing is larger (a^^ ~ 2 GeV). However, these errors remain 
substantial in the heavy quark region. We have implemented a proposal for reducing 
0{am) effects, which consists in matching the lattice quantities to their continuum 
counterparts in the free theory. We find that this approach still leaves appreciable, 
mass dependent, discretization effects. 
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1 Introduction 



One of the major goals of lattice QCD computations is the non-perturbative evaluation 
of hadronic Matrix Elements (ME's) which are relevant for strong and electroweak 
interactions. On the lattice, the ME's of bare operators are computed at finite UV 
cutoff (a is the lattice spacing) by simulating QCD at fixed bare coupling (3 = Q/qq. 
These ME's are functions of the finite cutoff. In order to extract physical information 
from them, the lattice operators require renormalization. For example, the p-meson 
decay constant /"-"^ is obtained from the ME (0| Vfc(0)|/3r.) as follows: 

2 

el-!L = Zv{mmPr) (1) 
Jp 

where nip is the meson mass and e)J is its polarization vector. T4 is a spatial compo- 
nent of the bare (lattice) vector current and Zy is its renormalization constant (RC). 
Zv connects the lattice theory at finite cutoff to the renormalized continuum theory. 
Analogously, the 7r-meson's decay constant /tt is determined from the ME of the axial 
current through: 

if^m^ = ZA(0|Ao(0)|7r) (2) 

where is the pion mass, Aq is the time component of the bare (lattice) axial current 
and Za is its RC. These are specific examples in which operator renormalization is an 
essential step for passing from lattice to continuum physics. This renormalization is 
an inevitable requirement in those cases where the ME's of bare operators diverge in 
the continuum limit (e.g. the ME's of the AI = 1/2 four-fermion operators). 
So far, three methods have been implemented to the computation of the Z's: 

• Lattice Perturbation Theory (FT) 

• Non-perturbative use of Ward Identities (WI's), either on hadronic states |11]-[20| 
or on quark states 21 1. 

• Non-perturbative (NF) renormalization on external quark and/or gluon states 
l2l-| 



Lattice FT can be applied to the calculation of finite and logarithmically divergent 
RC'sQ. The WI method is only applicable to finite RC's, whereas the NF one can be 
applied to all cases. The RC's, calculated to ©(c/g") in FT, suffer from 0(^0^"+ ^ 
systematic errors (in most cases n = 1). On the other hand, non-perturbative determi- 
nations of the RC's (WI or NF method), are subject to 0{a) systematic errors, since 
they are obtained from correlation functions computed at finite cutoff. Thus, physical 
quantities like those of eqs. (||) and (^), will suffer from 0{gfj^^^^^) and 0{a) systematic 
errors if the RC's are obtained in FT, and from 0{a) errors if they are computed non 



^PT cannot be applied to power divergent RC's, because of the possible presence of finite non-perturbative 
contributions (see, for example, pS], p8[ and references therein). 
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perturbatively. Therefore WI or NP estimates of the RC's are clearly preferable in 
principle. 

In recent years, there have been considerable efforts in refining the perturbative 
computation of the RC's. Boosted PT (BPT) has been extensively explored by Lep- 



age and Mackenzie [29|. It consists in judicious choices of the expansion parameter 
(motivated by the behaviour of the resulting perturbative series, by tadpole resumma- 
tion arguments and even by Mean Field Theory) which should bring the perturbative 
estimate of the RC's closer to their actual values. Even though the error remains 
0{gQ^^^^^), any of the several choices of BPT, proposed in [29|, might improve the 



asymptotic convergence of the perturbative series. 

On the non-perturbative side, Symanzik-type improved actions and operators have 
been constructed, in order to systematically reduce the finite cutoff error inherent in 
numerical computations. The implementation of the tree-level improved Clover action 
and operators reduces the discretization errors from 0{a) to 0{gQa) (more 

precisely, it removes all ©(a^g"" In" a) systematic errors, which are effectively 0{a) in 
the scaling limit, Qq ~ 1/lna). Recently, a non-perturbative extension of the tree- level 
Clover improvement programme has been proposed in ref. (3^, which corrects up to 
O(a^) errors. The general proposal of these works is directly applicable in the chiral 
limit (with the exception of the improvement of the vector current which has also been 
implemented at non-zero quark mass; see ref. |3^). Another approach to the non- 
perturbative improvement programme, which can also be applied with finite values of 
the quark masses, can be found in in refs. [34, 35]. 



Regarding the RC's, the main achievement of the tree-level improvement pro- 
gramme has been the reduction of discretization errors from about 20% to about 5% 
in computations with light quark masses [|l5[. The situation is more problematic with 
heavy quark masses, say am ^ 0.4. Even with a tree-level improved action, the dom- 
inant 0{gQam) error is appreciable. The most direct way of reducing these errors 
consists in reducing the lattice spacing (by going to larger /? values), but the increase 
in computational cost becomes prohibitive. Alternatively, one can try to reduce sys- 
tematically these errors at fixed cutoff. The proposal of refs. |3^] - [^] is such an effort. 
It consists in matching lattice quantities to their continuum counterparts at vanishing 
spatial momentum, ap, and large finite masses, am. The aim is to eventually remove all 
am corrections, leaving predominantly 0{ap) errors. This program has been explicitly 
carried out only in the free theory; we refer to it as the KLM approach. Speculations 
about higher order extensions can be found in [^ ]. 

This work is a thorough comparison of the various estimates (perturbative and non- 
perturbative) of the RC's of two-fermion lattice operators. We have tried to be ex- 
haustive in investigating the numerous improvement proposals found in the literature. 
In particular: 

• We have concentrated on the finite RC's of the vector and axial currents {Zy 
and Za) and the finite ratio Zs/Zp of the RC's of the scalar and pseudoscalar 
densities. 
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We have examined the influence of boosting to the 1-loop perturbative estimates 
of the above quantities. Several boosting recipes of ref. have been applied. 
The "optimal" scales of the boosted running coupling (the q* of ref. ||2^) has been 
obtained for Zy, Za and Zs/Zp. The q* values can be found in table |^. We find 
that all recipes of BPT give the same results within a few per cent (10% in the 
worst case). These results are summarized in table ^. Moreover, all BPT results 
differ from those obtained in SPT by up to 20%. When compared with the RC's 
obtained from WI's (see table BPT results are usually in better agreement 
than the SPT ones. 

The RC's Zy, Za and the ratio ZgjZp have subsequently been computed non- 
perturbatively, using the WI method, for several quark masses. Both the Wilson 
and the tree-level Clover actions have been implemented. Discretization errors at 
small quark mass, being 0{a) tend to be about 20% - 30% in the Wilson case. In 
the tree-level Clover case, being ©(^qq), these errors drop below 5%. However, for 
masses ma ^ 0.4, which are used in current lattice QCD simulations, they grow 
beyond 10%. By extrapolating our data in the quark mass, we obtain extremely 
consistent Clover results for Zy and Za at zero quark mass, shown in table ^. We 
have also checked that it is not possible to optimize these results by implementing 
different choices of improved operators, based on a different arbitrary admixture 



of the so-called Jf) and mo improving field rotations (see Subsection 5.2 for the 
definitions of these field rotations). 

In order to control discretization errors of increasing quark mass, we have applied 



the tree- level KLM matching prescription of refs. |36| - |^^. In the Wilson case, 
the large systematic error due to 0{am) effects is significantly reduced in some 
cases (see figure ||). There are cases, however, for which the KLM prescription 
fails to correct large 0{am) effects (see figure |2|). Moreover, a similar analysis, 
applied to the Clover case, revealed that KLM corrections fail to improve the 
situation. This is because the Clover correlation functions we have examined are 
characterized by large 0{gQam) errors in the region of heavy quark masses. Thus, 
the KLM corrections, being of O(a^), cannot absorb these discretization errors. 

The present work is intended as a thorough comparison of the various improvement 
schemes which have been widely used in present-day, state-of-the-art lattice QCD phe- 
nomenology studies. As such, it has not dealt with the more recent non-perturbative 
Clover improvement of |32, p3, 35 1, which appears to be the most promising way to 



reduce finite cutoff effects. 

The paper is organized as follows. In sect.|| we give some basic definitions. In sect, 
y we review the perturbative calculation of the RC's of two-fermion lattice operators. 



In subsec. 3.1 we discuss BPT and compare results for the RC's obtained from various 



boosting recipes. We also present the "optimal scales" q* of ||2^ for the operators of 
interest (some technical subtleties of this calculation can be found in the Appendix). 
In sect.Q we review the well known vector and axial WI's on the lattice, with particular 
emphasis on the conditions under which they can be used to obtain the RC's. We 
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have gathered in this section all WI's which are subsequently used in sects. |5| and 
where we discuss the practical implementation of the WI determination of RC's and 
present our results for Zy, Za and Zs/Zp. In sect. |^ we review the tree- level factors 
of the KLM-procedure for the Wilson correlation functions and derive their Clover 
counterparts. Finally, in sect. ^ we present our conclusion. 



2 Definitions and Generalities on RC's 

In this section we define the quantities of interest to this paper. Besides fixing our 
notation, we also review some general properties of the RC's of lattice operators, with 
particular attention to their dependence on the theory's parameters and renormaliza- 
tion scale. 

We consider bilinear quark operators of the form: 

0^{x)=i^{x)V^i,{x) (3) 

where ip{x) is the quark field, T denotes a generic Dirac matrix and are the genera- 
tors of the SU{Nf) fiavour group in the fundamental representation (a = 1, . . . , Nj — l). 
These generators satisfy: 



>J_ 

T'T 



where / represents the identity. Specific bilinear operators will be denoted accord- 
ing to their Lorentz group transformations: the scalar and pseudoscalar densities are 

denoted by = ilj{x)—ip{x) and P°'{x) = '4'{x)—"f^'il){x) respectively, whereas 

the local vector and axial vector currents by V^{x) = ipi^x) — 'y^'ip{x) and j4^(x) = 
A° 

'4'{x)~^lfil5'^ix)- Only such non-singlet bilinear quark operators will be considered in 
this paper. In the case of the vector and axial currents, we will also be using their 
point-split discretizations. 

We will deal with operator insertions in n-point Green's functions, with the operator 
usually inserted at the origin. For example the insertion of Or in the 2-point Green's 
function is 

Go{xi,X2) = (V(xi)O?(0)V;(x2)) (5) 
In momentum space the insertion becomes 

Go{pi,P2) = J dxidx2exp{ipiXi)exp{-ip2X2)Goixi,X2) (6) 
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and the corresponding amputated Green's function is given by 



^o{pi,P2)=S-\pi)Go{pi.P2)S-'{p2) (7) 

where S{xi — X2) = {ip{xi)ip{x2)) is the quark propagator and S{p) its momentum 
space counterpart. We will also be using Green's functions between external hadronic 
states. These are defined in an analogous fashion. For example, the 3-point correla- 
tion function {P{xi)V^{0)P^X2)) is a correlation function of three bilinear operators. 
Prom its Fourier transform in 3-dimensional space, and with the aid of the spectral 
representation formula at large euclidean time separations, the on-shell hadronic ME 
(7r(pi)|V^|7r(p2)) can be readily extracted. 

Everything defined so far is a bare quantity, which we assume to be defined in 
the lattice regularizationQ. We opt for the lattice regularization scheme proposed by 
Wilson which consists of a gluonic action Eol, 



S9 = -^T.Tr[Up + ul,] (8) 



and a fermionic part |41], 



X,fJ, 

+a4 5^V^(x)(Mo + -)V'(x) (9) 

where, in standard notation. Up is the Wilson plaquette, U^{x) is the lattice gauge link 
and go the bare coupling constant. The fermion field is a vector in flavour space. Its 
components will be denoted by u, d, s etc, or by ipi, ip2, fps, etc. according to notational 
convenience. The diagonal bare mass matrix is denoted by Mq and its elements by 
moi, ?n.02i nT-os etc. (or mou, mod, etc.). We use mo in the mass degenerate 
case. As is well known, for Wilson fermions, the quark mass, besides a multiplicative 
renormalization, is also subject to an additive one. Defining m = mg — rric, the 
renormalized mass will be given by: 

ruR = Zmin = Zm[mo - mc] (10) 

with Zm the multiplicative RC. The chiral limit is then m — > 0; i.e. rriQ — > mc- At 
tree-level, nic = — 4/a. For the diagonal elements of the matrix Mq — mc we will use 
either mi, m2, ms etc., or m„, m^, m^ etc., according to notational convenience. 

All operators studied in this work are subject to multiplicative renormalization; 
cases involving mixing of equal- or lower-dimensional operators (e.g. the four-fermion 



^This means that the integrals of eq. (g) are really sums (a* X^ici which run over all lattice sites, 
labelled by xi, X2, etc. The use of integrals instead of sums and Dirac functions instead of Kronecker 
symbols is a question of notational convenience, which should be transparent. 
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operators driving the AS* = 2 or AI = 1/2 transitions) will not be considered. Thus, 
the renormalized operator is given by 



(11) 



where gR is the renormalized gauge coupling and /x the renormalization scale. Note that 
in the bare quantities, the bare mass mo has been traded off for the more convenient 
(from the point of view of the chiral limit) subtracted mass ra. A specific, if arbitrary, 
choice of renormalization scheme and scale are implied. Suppressing the dependence on 
couplings, masses, scales and cutoffs for notational simplicity, we list other relationships 
of interest between bare and renormalized quantities (denoted in general by "hats"): 
the renormalized quark propagator is 



Let us now examine the dependence of the RC's Zq on the bare parameters of the 
action. Multiplicative Zq^s are dimensionless quantities. Consequently, their functional 
dependence on the bare parameters and the cutoff can in principle be only of the form 
Zo{gQ:0,fn^o,lj)- For the specific case of finite RC's (e.g. those of the vector and axial 
currents) there is no dependence on the renormalization scale /z. In other words, the 
anomalous dimension is zero; this will be discussed in detail in the following section. 
Therefore, we are left with a possible dependence of the form Zo{gQ,am). However, 
there can be no dependence of Zq on am. Regular terms such as positive powers of 
am, which may appear in the calculation of the bare Green's function Ao, will drop 
out of Zq, since RC's are strictly defined in the continuum limit a — > 0; c.f. eq. (p^. 
This leaves us with singular terms, such as In(am) and negative powers of am. But 
since we require the existence of the chiral limit (m 0) in QCD, such a dependence 
must also be ruled out. This exhausts all possibilities, and renders Zq independent 
of am. In conclusion, the finite RC's Zy, Za and the ratio Zs/Zp, examined in this 
paper, are of the form Zo((?o)-0 

■^Dimensionless RC's which are not finite have the general dependence Zo{gQ, am, a/i). Provided that the 
renormalization scale /x is chosen to be greater then the quark masses (/i 3> m), there can be no dependence 
on am (same argument as above). The requirement of a well defined chiral limit allows only regular terms 
in m/ fi, which however have been assumed negUgible. Thus, in this case the RC's have the form Zo(g§, a/x) 



S{p) = lim[Z^S{p)] 



(12) 



where Z^ is the quark field RC. From the above definitions of Zq and Z^, with the 
aid of eqs. (|5|) - (0), we find for the renormalized Green's functions: 




60(^1,^2) = lim[Z^ZoGo{pi,P2)] 
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3 Evaluation of RC's from Perturbation Theory 



In this section we review the basic concepts of the perturbative evaluation of the RC's of 



lattice operators and examine the proposals of ref. |29| which aims at the improvement 
of the convergence of the perturbative series in lattice PT. 

The standard perturbative determination of the RC Zq of a bilinear quark operator 
Or (flavour indices are suppressed, unless necessary) requires the calculation of the 
quark self- energy and the amputated Green function Ao(j)) of eq. (|^) between external 
quark states, both at momentum (see, for example, [^). Here we will work with the 
projected quark propagator Vs(j>)-, and projected amputated Green function Vo(j>)-, 
defined as 



dp,. 



(14) 



where the trace is over spin and colour indices and Pq is the Dirac matrix which 
renders the tree-level value of To{p) equal to unity (i.e. it projects out the nominal 
Dirac structure of the Green function Aq (p) ) : 



Ps 
Pv 



r7M 



Pp = 75 

Pa = 



Prom eq. (|T^ we see that rs(p) renormalizes like 



The renormalization of Po(p) is identical to that of Ao(p); c.f. eq(|l3| 



fo(p) = lim \z^^ZoAoip) 



(15) 



(16) 



(17) 



In the last two equations, the limit corresponds to the removal of the UV cutoff (e.g. 
a on the lattice). The choice of Ps(p) and To{p), being scalar quantities, simplifies 
the discussion. Projected Green's functions have been implemented in recent works 
(e.g. ref. p^). Another simplification consists in working with zero quark mass. This 
is justified, as we have shown in the previous section, by the lack of mass dependence 
of Zq, in the limit m. 

Upon passing from formal definitions of renormalized quantities (like eq. (0)) to 
their specific evaluation, we need to implement a regularization. The standard practice 
in continuum calculations consists in choosing a Dimensional Regularization {DR) 
scheme, such as Naive Dimensional Regularization (NDR), 't Hooft - Veltman (HV) 
etc. The expression for To, evaluated at 1-loop PT, can be written in the general form: 



rg^(p,<?o,e) 



1 + 



(47r)2 



1 



(7r- - 7rln(p//i) + Cj 



+ Oie) 



(18) 
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where 7r is an anomalous dimension and C^^ is a finite constant, which depends 
on the regularization scheme and the gauge chosen. The factor 1/e in eq. ( p!^ ) is an 
abbreviation for 

\ = ^+ ln(47r) - jE (19) 

where e = (4 — D)/2 and je stands for Euler's constant. In DR we work in D 
dimensions, where the original bare couphng go has dimension e. Here the scale fi is 
introduced to render the bare couphng goifJ-) dimensionless. The /U-dependence of the 
r.h.s. of eq. ([l^ is apparent. 

In order to renormalize Tq, we are in principle free to impose any renormalization 
scheme. The MS scheme, at 1-loop, amounts to removing the 1/e divergence. Since 
the RC of the projected Green function Tq is given by Zr = Z^^Zq (see eq. ([T7|)), 
this implies for the value 



,MS,DR 



(5o(At),e) 



1 



1 

(4^)2 ^^e 



(20) 



Consequently, the renormalized Green function is given by 



1 + 



(4^)2 



-7rln(p//i)^ 



C7^^) 



(21) 



where, to this order, we are free to replace go by 5^(^)5 the MS renormalized coupling 
constant. 

The same calculation can be repeated on the lattice. Now the UV cutoff is provided 
by the inverse finite lattice spacing and thus the 1-loop calculation yields: 



r 



LAT 
O 



{p,go{a),a) 



1 + 



go (a J 
(47r)' 



-(-7rln(pa)^ + Cf^^) 



+ 0{a) 



(22) 



where go is the bare coupling of the lattice action. Note that at 1-loop the anomalous 
dimension 7r is regularization independent. 

The renormalization scheme can again be chosen at will; often the MS is chosen 
also on the lattice. This seemingly unnatural choice (the MS is closely linked to con- 
tinuum DR) has a few advantages. For example, ME's of effective Hamiltonians, once 
calculated non-perturbatively on the lattice, must be renormalized and combined with 
perturbatively calculated Wilson coefficients, in order to obtain physical amplitudes. 
The renormalization-group invariance of these amplitudes is guaranteed only if the 
Wilson coefficients and the RC's are calculated in the same renormalization scheme. 
Since the former are often known in the MS scheme, this scheme is also preferred for 
the calculation of the latter. This choice then corresponds to the requirement 



lim 



Z 



MS, LAT 



{l^a,go{a))To^^{p,go{a),a) 



(23) 
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where again the lattice couphng go{a) should be traded for the MS renormalized cou- 
pling constant g^{fi). This point of principle is of limited relevance for a 1-loop 
calculation. From eqs. (H), (H) and (H), the following RC is obtained: 



Zp {fia,go{a)) = 1 + 



(47r)2 



7rln(/xa)^ + C^^-C^ 



^LAT 



(24) 



We now recall that the RC of the amputated vertex Tq is = Z^^Zq- The RC of 
the quark field Z^ can be calculated, from the renormalization of the quark propagator 
S{p), with an analogous procedure (c.f. eq. (|l^)). The result is 



{^la, go{a)) = 1 + j^-^ 
Combining eqs. (24) and (|25|) we obtain 



iLAT 



yMS,LAT 



(/xa,go(a)) = 1 + 



50 (a 



i2 . 



(47r)^ 



(25) 



(26) 



where 



70 
Ar 
As 



7r + 7s 
r<DR r<LAT 

Op — Lyp 

r'LAT 



(27) 



It is this RC (with this choice of renormalization condition) which is usually denoted 

by Zq in lattice PT calculations. The dependence of z^^'^^'^ on the coefficients C^^ 
and comes from the choice of the MS renormalization condition (see eqs. (^Tj) 

and (|2^)) whereas its dependence on Cp"^"^ and 0^/^'^ from the lattice regularization 
(see eq. (p^)). Two perturbative calculations are thus necessary, one in the continuum 
for the C^^'s and one on the lattice for C^'^^^'s. The presence of the C^^'s on the 



r.h.s. of eq. (24) is sometimes called regularization dependence of the renormalization 
scheme The reader is referred to ref. Q for the explicit calculation of the RC's 

of several bilinear operators with the Wilson action. The corresponding results for the 
Clover action appeared in ^. 

The PT method is applicable to finite and logarithmically divergent operators, the 
renormalization of which does not require a subtraction of lower dimensional operators. 
Model examples of operators with logarithmically divergent multiplicative RC's are the 
scalar density S{x) = ^i{x)ip2{x) and the pseudoscalar density P{x) = ipi{x)j5tp2ix). 
These are the cases best suited to our general discussion so far. 

The RC's of the vector current V^{x) and the axial current ^^(2;) are also calcula- 
ble in PT on similar lines. Being (partially) conserved quantities however, they have 
some peculiarities. Their (partial) conservation guarantees zero anomalous dimension 
(i.e 7y = and = 0). This result is regularization independent (e.g. valid both in 
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the continuum (DR) and the lattice (LAT). It imphes that the current renormaUza- 
tion constants are finite, scale independent quantities. Moreover, the following useful 
properties must be taken into consideration: 

• The vector current conservation survives the continuum DR. Thus it satisfies a 
WI which implies that z^^'^^ = 1. This is an exact result, valid at all orders of 
PT. Therefore, at 1-loop V^{x) has zero finite constant; i.e. C{^^ + = 0. 

• As we will see in the next section, the local lattice vector current is not a conserved 
quantity. It has a non- vanishing 1-loop contribution (i.e. Cy^^ + C^'^'^ 7^ 0)- In 

other words, Z^^'^^"^ (qq) ^ 1. 

• We will also show in the next section that on the lattice one can define an ex- 
tended conserved vector current, V^{x), which is also trivially renormalized with 
^MS, LAT _ ^ Thus, also for this current we have C^^^ + C^^^ = 0. 

• For the axial current the argument is almost the same. However, some further pe- 
culiarities arise due to the breaking of the chiral symmetry by the regularization. 
In the continuum, the finite 1-loop contribution to the renormalization constant 
may or may not vanish, depending on the detailed choice of the regularization 
scheme. Thus, for NDR we have C^^^ + C^^^ = 0, implying z^s,iVDi?, ^ 
On the other hand, with HV for example, we have C7f ^ + C^^ / 0, implying 
^MS,//y^^2-j _^ ^ present work, we will always imply NDR when referring 
to DR. 

• We will show in the next section that the lattice regularization breaks chiral 
symmetry. Thus on the lattice C^"^"^ -|- C^^^ ^ (for any definition of the bare 

lattice current, local or extended). In other words, Z^^'^"^"^ {g^) ^ 1 

We are also interested in the finite ratio Zs/Zp. Its finiteness relies on WI's (see 
next section) and implies that 75 = 7p. With the specific choice of NDR as the 
continuum dimensional regularization, we also have Cg^^ = Cp^^. 

3.1 Boosted Improvement of the PT Estimates of RC's 

The perturbative calculation outlined above is only approximate, because of the trun- 
cation of the perturbative series. At present, most lattice RC's are only known in 
PT at 1-loop. Thus, they suffer from O(5o) systematic errors. These errors must be 
born in mind when comparing the RC's, calculated in PT, to their non-perturbative 
estimates. In this subsection we discuss these systematic errors and comment on the 



proposals for improving the convergence of the perturbative series of ref. [29|. These 
proposals are known as Boosted Perturbation Theory (BPT). 

We have already pointed out that, in principle, at 1-loop there is an ambiguity in 
the choice of couplings {go{a), go{fJ-), 5ms (^)) ^° used in eq. (^). This ambiguity 
is lifted in (and up to) 2-loop PT. At 1-loop however, a judicious choice of expansion 
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parameter may be important in practice. The lattice couphng g^ia) has been shown to 
be a bad expansion parameter in ref. |29|, where several improved (boosted) couplings 



have been proposed, and several observables have been calculated in PT with these 
couplings. 

In the specific case of Zq, the Standard Perturbation Theory (SPT) result is given 
as a series in the bare (lattice) coupling go 

Zo{ai^,g^o{a)) = 1 + C7i(a/i)g2(a) + C2(a/x)g^(a) + . . . (28) 



To first order, the above equation reduces to eq. (pq), with Ci(o/i) = [70 ln(a;u) + 
Ar + As]/(47r)^. According to (2^, perturbative series have a better behaviour if 
expressed in terms of a renormalized coupling gy{q*) at a suitable scale q* (to be 
specified shortly): 

Zo{afi,gUql) = I + Ci{api)gl{q*) + C2{a^,,aq*)9Uq*) + ■ • • (29) 

The claim is that the perturbative estimate obtained from the boosted series, truncated 
at low orders, is closer to the non-perturbative result than its SPT counterpart. 



The renormalized coupling gv{q) in eq. (p9) can be chosen at will. In ref. |29] it is 



suggested to define it from the heavy quark potential V{q) at momentum transfer q: 

V(,) = -'4^ (30) 



Moreover, according to |29| and |44], the scale q*, which is most appropriate to the 



expansion (g9|) of Zp, is fixed by: 

where f{q, afj,) is the function entering in the loop integral which defines the constant 
Ci{afi) of eq. ([29|): Ci(a/x) = / d'^qf{q, a/^).0 From the above expression we readily see 
that in general q* depends on the renormalization scale ^u. Here, however, we are only 
interested in the finite quantities Zy, Za and Zs/Zp, for which the Cj's of eq. (28) are 
/i-independent. Consequently, also q* does not depend on //. 

The connection between the bare coupling gQ and the renormalized coupling gy 
defined in eq. ( |30|) is obtained by calculating the heavy quark potential in lattice PT. 
At lowest order one finds 



9viQ) = 5o(«)[l - gi{a)2l3oHaq/c) + Oig^{a))] (32) 

"^Here we are dealing with multiplicative renormalization, for which the above recipe is fairly non ambigu- 
ous. However, if the operator mixes with others under renormalization, the recipe is ambiguous, since one 
can either evaluate all mixing constants at the same scale q* or obtain each mixing constant at a separate 
scale. 
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where Pq = ll/(47r)^ and c is given by: 



There are several possibilities, suggested in for extracting gy{q*) from eq. (p2D. 
They can be classified as follows: 



1. Purely Perturbative Boosting [PPB]: This method consists in using eq. (|3^ ) in 
order to identify the couplings at scale q = c/a: 

gl{c/a)=gl{a) + 0{9',) (34) 

Then one can perform 2-loop Renormalization Group (RG) running of gY{c/a) down 
to the scale q* in the standard fashion: 

1 - 1 +2;%ln(5?:)+ftl„f4<fM) (35) 



9viQ*) 9v{c/a) c /3o V gyiq* 



with /?! = 102/(47r)^. Thus, having calculated q* (from eq. (|3l|)) and gy{q*), we derive 
Zo from eq. (|29|). We call this estimate PPB. 

2. Monte Carlo - Perturbative Boosting [MC-PB(ln □)]: Again this method is based 



on the relationship between gy and g^ of eq. (|32|) . But now the boosted coupling con- 
stant is extracted from the non-perturbative (Monte Carlo) value of a short distance 
quantity, such as the average plaquette |4|]. The logarithm of the plaquette, known 



as a perturbative series in g^, is expressed, with the aid of eq. (32), as a series in gy. 



The leading coefficient of this series is an integral of the kind / d qfn{q), which gives 
rise, through the criterion of eq. (|3l[), to the scale (?□, appropriate for this expansion. 
Prom ref. we have q^ = 3.41/a and the expansion: 

H^TriUu)) = -lgUQn)[l - 9.46 • lO'^gUQu} + Oig^ (36) 

The series of the plaquette's logarithm has been preferred to that of the plaquette itself 
because, as claimed in ref. [p9| , its perturbative expansion is better behaving. Using 
the Monte Carlo result for the average plaquette, eq. ( p6D is solved for 5y(g'n). 2-loops 
RG running subsequently yields the desired gy{q*)- With this value for gy{q*), we 
derive Zq from eq. (|2|). We call this method MC-PB(lnn). 

3. Monte Carlo- Perturbative Boosting [MC-PB(n)]: Yet another proposal of p9| 



consists in combining eqs. (32), at the scale q = it /a, and the perturbative expansion 
of the plaquette 

^TriUn) = 1 - IgU^/a) + 8.71 • W-^gU^r/a) + ©(g^-) (37) 



to obtain: 



gU^/a)[l - 0.513^5^(7r/a) + 0{gt.)] = (38) 
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Wilson 


Clover 


aq*{Zv) 


aq*{ZA) 


aq*{Zs/Zp) 


aq*{Zy) 


aq*{ZA) 


aq*{Zs/Zp) 


2.4 


2.6 


1.9 


2.7 


1.2 


3.2 



Table 1: The values of the scale aq* , characteristic of the boosted perturbative evaluation of 
Zy, Za and Zs/Zp. The scales are shown for both the Wilson and Clover actions. 



Once (7y(vr/a) is obtained from the above expression, two loop RG running yields 
gy{q*) and subsequently Zq- We call this estimate MC-PB(n). 

4. Naive Boosted Perturbation Theory [NBPT]: A particularly simple choice of boo- 
sted coupling 1 45] is = g^/ {^Tr{Un)) and the subsequent substitution of by g^ 
in eq. (24). This recipe (a simplification of eq. (|38[)) is denoted by NBPT. It has been 
implemented, for example, in refs. |46| and ||4^. As argued in ref. Mean Field 

arguments can also be used in order to support this prescription. 

Having classified the various boostings, we now comment on them. All these meth- 
ods should be in principle equivalent to lowest order in PT. From the field theoretic 
point of view, PPB is a perfectly legitimate operation, as it amounts to a different 
choice of expansion parameter. The two MC-PB prescriptions may be even superior 
to PPB in practice, since they possibly incorporate some non-perturbative and higher 
order effects. In order to implement these recipes to the perturbative evaluation of the 
RC's of interest, we have calculated the appropriate scale q*; this is different for each 
RC. The values of q* for Zy, Za and Zs/Zp are collected in table [l|. Some technical 
details concerning their evaluation are gathered in the Appendix. 

In table |2|, we collect the values of Zo, obtained from the various boosting recipes 
listed above. We see that the variation of the boosted coupling, due to the choice of 
boosting recipe, does not exceed 20%, and that of the various RC's is about 10%. These 
variations reflect the systematic error introduced by the truncation of the perturbative 
series. In the next section, we will compare boosted results to the ones obtained 
non perturbatively from WI's. For this comparison we will use the MC-PB(lnn) 
prescription which is considered the optimal choice by the authors of ref. |29]. As will 
be shown in table ^, such a comparison reveals a general trend for BPT to reduce the 
magnitude of the ©(^g) systematic error. However, in some cases this error remains 
appreciable (more than 10%). 



4 Lattice WI's with Wilson Fermions 

In this section we review the WI's which follow from the chiral flavour symmetry of 
the QCD action. They allow to compute the finite RC's of lattice operators in a 
non-perturbative way. Several properties of these RC's can be derived from the WI's. 



Although the results presented here are well known (see [48, ^), we give particular 
emphasis to the conditions (e.g. continuum limit, chiral limit etc.) under which these 
results are exact. This is essential to the understanding of the sources of systematic 
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Wilson 


- /3=6.0 










Zv 


9' Za 


9' 


Zs/Zp 


SPT 


1.00 


0.83 


1.00 0.87 


1.00 


1.10 


PPB 


1.81 


0.69 


1.77 0.76 


1.93 


1.25 


MC-PB(lnn) 


2.14 


0.63 


2.08 0.72 


2.32 


1.34 


MC-PB(n) 


1.97 


0.66 


1.92 0.74 


2.13 


1.29 


NBPT 


1.68 


0.71 


1.68 0.77 


1.68 


1.20 






Clover 


- /3=6.0 








a' 


Zv 


Za 


9' 


Zs/Zp 


SPT 


1.00 


0.90 


1.00 0.98 


1.00 


1.20 


PPB 


1.75 


0.83 


2.25 0.96 


1.67 


1.42 


MC-PB(lnn) 


2.05 


0.80 


2.81 0.95 


1.95 


1.54 


MC-PB(n) 


1.90 


0.81 


2.52 0.96 


1.81 


1.48 


NBPT 


1.68 


0.83 


1.68 0.97 


1.68 


1.43 



Table 2: Wilson and Clover action RC's obtained with different boosting procedures (see text 
for notation). The RC's are shown at (3 = 6.0. 



error in practical applications of the WI's, which form the basis of our analysis. The 
reader familiar with refs. p8| , |4^ ) may skip this section. 

We consider local SUl{N j)xSUr{N f) chiral transformations of the fermionic fields: 



5il}{x) 
6'tlj{x) 



A" 



avUc 



(39) 



and examine the WI's derived separately in the vector and axial cases. 



4.1 Vector WI's on the Lattice 

With degenerate quark masses (i.e. Mq in eq. (^ proportional to the unit matrix) 
global vector transformations are a symmetry of the action. From the corresponding 
local transformations (eqs. ( |39[ ) with = 0) the following vector WI can be derived 
HI: 

.^ SO{xu...,Xn) ^ = a^EV^(0(xi,...,rrjf^(x)) + 
oay[x) ^ 

+a4(0(xi, . . . ,xjVS(x)[y , Mo]V'(x)) (40) 

where 0{xi, . . . , Xn) is any multi-local operator consisting of quark and gluon fields at 
different space-time points (xi / X2 / ...x^). aV^/(x) = (/(x) — /(x — /x)) is an 
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asymmetric lattice derivative, and 



V^{x) = \[i>{x){^,-l)U,{x)^^l,{x + iJi) + 

+^(x + ^)(7^ + l)C/t(x)yV(x)] (41) 

is a point-split vector current. By keeping the point x separate from the points 
xi, . . . , and by taking the limit of degenerate bare quark masses, we see that eq. ( ^0[ ) 
leads to the conservation of the point split lattice vector current, V^V^(x) = 0. This 

also implies that the standard local vector current VJ^(x) = ■ip{x)—^f^ip{x) is not con- 
served on the lattice. 

For the specific choice 0(xi,X2) = '^p{xl)^p{x2), the vector WI (4C) becomes 

Y,V^{^l^{xl)V^{x)i;{x2)) = {ij{xMx)[Mo, ^Mx)i^{x2)) + 

+6{x2 - x){'ip{xi)ij{x2))Y - ^i^i - Y(V'(a:i)V^(a;2)) (42) 

Let us consider the case of degenerate quark masses (Mq oc 1) and place the vector 
current at the origin (x = 0) . By Fourier - transforming the above equation and ampu- 
tating the corresponding Green's function = {4'{xi)V'^'°'{0)'ip{x2)) (c.f. eqs. (|5|), 
(^) and (0)), we obtain 



E 



1 - exp{-iaq^) 



Ar(f + |.P-|) =-4 + 1)" ^4^11-^1 (43) 



This is the lattice version of the vector WI. In the continuum limit, it is a formal 
expression relating bare quantities which are divergent. In order to render them mean- 
ingful, it is essential to go over to renormalized expressions. With the RC's strictly 
defined in the limit of vanishing cutoff, we express the bare quantities of eq. ( |4^ ) in 
terms of their renormalized counterparts, given in eqs. (H) and (111). Thus, in the 
limit a ^ 0, we obtain: 

E>*.^,'Ar (p + |,P-f) =-<P + f)"f + f<p-|)" (44) 

However, the definitions of the RC's of eqs. (|l2|) and (|^), used to derive the above 
result, are only formal. A concrete definition implies the choice of a renormalization 
scheme. In the case of the vector current, it is crucial to impose that the renormalized 
operator satisfy the nominal WI of the continuum theory. From this requirement 
on eq. (p^ ) we obtain that the RC of the point split conserved vector current (|4l| ) is 

Hiii: 

Z~ = 1 (45) 
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This result, which is exact in the limit a ^ 0, guarantees a proper definition of the 
vector charge and the validity of current algebra. Thus, the vector symmetry of the 
theory survives renormalization. 

The local vector current V^(2;) = ^(x)— 7^'i/'(^) is not conserved on the lattice. 
However, it can be shown that its ME's differ from those of the conserved current by 
finite contributions. The argument goes as follows: first we express the the conserved 
current as 

y^{x) = + |[Vi(x)(7M - 1) Y V'(^) + ^(x) (7m + 1) yV'(x)] (46) 

where we have used the lattice asymmetric covariant derivatives 

a Vfi ip{x) = Uf,{x)ip{x + /i) - Tpix) 

aij{x)Vf^=i^{x + fi)Ul{x)-ij{x) (47) 

The second term on the r.h.s. of eq. ( |4^ ) is a dimension-4 operator (we call it A^) 
multiplied by the lattice spacing. For definitiveness, we now consider the Green's 
function A~(p) of eq. (0). From eq. (^) it follows that 

A~{p) = Av{p)+aAA{p) (48) 

The term aA/\{p) vanishes at tree-level in the naive continuum limit. Beyond tree-level, 
however, this term contributes, due to power divergences induced by the mixing with 
lower dimensional operators. The mixing with operators of equal dimension gives at 
most logarithmic divergences In(ap); thus in the continuum limit such contributions 
vanish like aln{ap). As far as lower dimensional operators are concerned, AJ^ only 
mixes with the dimension-3 vector current V^. This gives rise to power divergences a~^, 
without logarithms^. Thus, the term cAa gives finite contributions, which combine 
with the renormalization of Ay to give the known non-renormalization exact result 
Z~ = 1 of the l.h.s. of eq. (^). Consequently, the term Ay on the r.h.s. has a finite 
renormalization: 

Zvigi) + 1 (49) 

Being finite, Zy can only depend on the coupling g^. Note that eq. (and any other 
vector WI) can now be expressed in terms of the non-conserved local current; we simply 
substitute Ay by ZyAy. This is the basis of the non-perturbative WI calculation of 
Zy, detailed in sect. |5[ 

One can also construct useful WFs based on hadronic correlation functions. For 
definitiveness, we consider the case in which the operator 0(xi, . . . ,x„) of eq. (|40| ) is 

^The absence of logarithms in such mixing has been explicitly demonstrated in ref. |]50| for the axial 
current at all orders in PT; the situation is analogous for the vector current. 
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0{xi,X2) = P"^^(xi)P^^(x2), where P^^ = V'i75^2- With x ^ xi,X2, so that contact 
terms do not contribute, eq. (^) becomes 

Y^V%P^\x{)Vl\x)P''\x2)) = {m2-m3){P'\xi)S^'{x)P^\x2)) (50) 



In analogy to eq. (44), we can use eqs. (10) and (|TT]) in order to express the above 
identity in terms of the renormahzed densities, current and masses. By requiring that 
the renormahzed quantities obey the nominal vector Wl, we then find that the product 
of the quark mass times the scalar density is renormalization group invariant. Thus, 
we obtain for the RC of the scalar density S{x): 



(51) 



4.2 Axial WI's on the Lattice 



Far less straightforward is the implementation of the axial symmetry with Wilson 
fermions, because of the presence of the chiral symmetry breaking Wilson term in the 
action. This topic has been dealt with in refs. |48, 49|. The basic idea is that, by 
imposing suitable renormalization conditions, PCAC is recovered in the continuum. 
The axial WI, obtained from eq. 

,50{xi,...,Xn). 4 



with ay = 0, is 



5a\{x) 



Y^V'i{0{x^,...,Xn)Al{x)) 



-A" 



-a\0{x^, . . . , x„)Vi(x){y , Mo}75V'(^)) - a\0{xi, x„)X'^(x)) (52) 



where A'^(x) is a bilinear point-split axial current given by 

^m(^) = ^[V5(^)(7A.75)f/M(x)y^/'(x + /i) + 

A" 

+V(a; + f^){lf.l5)UUx)—ip{x)] 



(53) 



The term X"" in the above Wl is the variation of the Wilson term under axial trans- 
formations: 



where 



X^ix) 



1 

--a 
2 



V'(a;) V il^{x) + ^(x) V —lbi^{x) 



(54) 



a? V V'(x) = \u^{x)ip{x + + Ulix - fi)ip{x - /x) - 2V'(x) 
a^i^{x) V = [i){x + p)Ul{x) + ^(x - n)Uf,{x - - 2^(x) 
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(55) 



Unlike the vector current case, X° cannot be cast in the form of a four-divergence. 

X" is a dimension-4 operator which, in the naive continuum hmit vanishes, being of 
the form (a x dimension-5 operator). However, it has divergent ME's beyond tree-level. 
Its mixing with operators of equal and lower dimensions, worked out in ref. |4£], can 
be expressed as follows: 

X'^ix) = X-{x) + i;{x){^,Mh,i;{x) + {Zj- l)V^iA;{x) (56) 

where naive dimensional arguments tell us that the mixing constant Z^{gQ , am) is 
finite, whereas M{gQ, Mq) diverges linearly like a~^. Logarithmic divergences have 
been shown to be absent at all orders in PT ||5^. The mixing constants Z^{gQ,am) and 
M((7q,Mo), and therefore X"(x), are defined by requiring that the renormalized axial 
current satisfies the nominal (continuum) axial WI. This amounts to a renormalization 
condition which ensures that the renormalized theory has the proper chiral symmetry. 
For example, with 0{xi,X2) = ^/'(xi)V5(x2), eq. (|5^) becomes 

V^A'^{x,)Z^A;{x)i;{x2)) = 

= {i,{xMx){^, (Mo - M)}j5i^{x)i^{x2)) + (V(xi)X"(x)V;(x2)) - 

- A'^ A*^ - 

-6{x - X2)(V'(a;i)V'(a;2)75) Y - - a^i) Y(75V'(a^i)V'(a;2)) (57) 

where X"- has been replaced by X"" . We now take the continuum limit of the above 
equation, with the aim of obtaining from it the nominal axial WI for the renormalized 
quantities. We immediately see that this identity can be recovered from eq. ( p7| ) 
provided that the following renormalization condition is imposed on X : 

lim{iJ{xl)X''{x)^jJ{x2)) =0 (58) 

Moreover, the chiral limit is to be defined as the value Mc of Mq for which the difference 
Mc — M{gQ, Mc) vanishes. Thus, in the chiral limit, eq. ( |57| ) becomes 

V^(V(xi)ZjI^(:r)VS(x2)) = (59) 

- A"' A*^ - 

-5{x - X2){'4>{xi)i){x2)^^)— - 5{x - xi) — {"f^il:{xi)i}{x2)) 

Expressing the quark fields in terms of their renormalized counterparts, we immediately 
see that the nominal axial WI is recovered up to vanishing cutoff effects, provided that 
we interpret the product ^4^ = Z^A^j^ as the renormalized axial current. Note that 
Zj((7o , a'Ti) , having been fixed in the continuum limit by eq. (|58|), does not depend on 
am (only a regular dependence on am permits a well defined chiral limit). Thus, it is 
a finite renormalization of the form Z^lg^). 

The above arguments rest on the assumption that the chiral point Mq, defined 
after eq. (|58|), is consistent with other definitions (e.g. the vanishing of the quark 
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mass in PT or of the pion mass). This has been verified both in 1-loop PT |49] and 
non-perturbatively (e.g. ]l6|). 

In practice it turns out to be more convenient to work with the lattice local axial 
current We can show, in a fashion analogous to the case of the vector current 

(c.f. the power counting argument based on eqs. (^) - (^)), that A'^{x) has a finite 
normalization constant Za- Thus we have A'^ = liina—>o[Z ^A*^] = lima^o[-^A^Ji]- From 
now on, the combination Z^A'^^ will always be substituted by ZaA'^ wherever it appears 
in a WI. Also analogous to the vector current case is the lack of mass dependence of 
these constants. We therefore have: 

Z~{gl) , ZA{gl) + 1 (60) 

Hadronic axial WI's are also very useful. In eq. (^) we consider the operator 
0{x\) = P'^^{xi) = '>p2{xi)^5il^i{xi) and x / xi, so that contact terms do not con- 
tribute. The resulting WI is 

ZAj2^^Ml'ix)P'\x^)) = {x'\x)P'\x,)) + 

+ [moi + mo2 - Ml -m2](P^^(x)p2i(a;^)) (61) 

where mi((7o, Mq) the i*^ diagonal element of M((7o, Mq). Note that we have replaced 
X by X according to eq. (|56|). The condition (|5^ ) implies that correlation function 
of the operator X vanishes in the continuum limit, except for localized contact terms. 
Thus, when the points xi and x are kept separate, {X{x)P{xi)) vanishes for a ^ 0. 



We now write in eq. (61) the renormalized quark mass as: 

dfn 



niR = Zm[mo - m{mo)] = Zm[mo - mc 



dmn 



imo-mc) + ...] (62) 



mc 



where we have used the property rn{mc) = mc- Upon expressing eq. (|6l]) in terms of 
renormalized quantities and requesting the validity of the axial WI in the continuum 
limit, we obtain for the RC Zp of the pseudoscalar density: 

Zp = l/Z^ (63) 



Note that from the WI of eq. (61) we can only determine the combination 



12-1 _ _ EuV?^(Al2(x)p21(2;^)) 

2p = Z^ [moi + mo2 - mi - ma] = ^pi2(^)p2i(^^)) (64) 

but not Za- Again, the above is true up to vanishing cutoff effects. 
A very useful WI can be derived taking 0{xi^X2) = ^|^(xi)VT(x2): 



ZAV^AKi^)Al{x,)v^{x2)) 

^(x){^A",(Mo-M 
+ir^H{x - x^){V,\xi)V^{x2)) + irH{x - X2){Al{xM%X2)) (65) 



{i>{x){\\\ (Mo - M)}^^^{x)Ai{xi)V;;{x2)) + {X\x)Ai[xi)V;;[x2)) + 
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where again use of eq. (^) has been made. By symmetry arguments the contact terms 
to which X gives rise must have the form 

(X''{x)At{xi)V^'{x2)) = -icir'H{x - xi)ly,\x,)V^{x2)) + 

+ic2r''5{x - X2){Al{xi)A%^2)) + ... (66) 

where the dots stand for locahzed Schwinger terms, which will vanish after integration 
over X (we always keep xi / X2). We now proceed as follows: (I) rewrite eq. (|65[), 
using eq. ( |66|) and expressing all bare quantities in terms of the renormalized ones; (II) 
recall that the RC Z„ of Mq - M is the inverse of the RC Zp (c.f. eq. (H)); (III) 
require, as always, that the renormalized quantities obey the nominal axial WI. Thus 
we obtain 

C2 = ^-l (67) 
The above conditions ensure the recovery of the canonical WI in the continuum limit. 



This WI is the basis of the non-perturbative calculation of Za in the fashion of refs. |11 



17|, where eqs. (^) - (67) are combined into 

{[V'iAl{x)-i^{x){^,p]^,^P{x)]Al{x,)V^{x2)) = 

+ir'H{x - xi)(F/(xi)F;(x2))|| + ir'^Six - x2){Al{xM%^2))^ (68) 

with p = Z^^{Mq — M) a matrix generalization of eq. (|6^. Performing the integration 
over X, and integrating over xi (recall that xi ^ X2 is necessary in order to eliminate 
Schwinger terms) we arrive at: 

dx J dxMx){h'',p}j5i'{x)Ai{x,)Vp'{x2)) = 

-i^r"'/ dx^{V,\x^)V^{x2))-i^r'' j dxMi{xi)A%X2)) (69) 

Note that this equation is only valid away from the chiral limit, where the integral over 
X of the total divergence V^A^(x) vanishes. At zero quark mass, the term containing 
the total divergence contributes because of the presence of massless Goldstone bosons. 

The same arguments may be repeated for the operator 0{xi,X2) = S^{xi)P^{x2)- 
By imposing that the renormalized operators S*-^ (xi) and (X2) satisfy the (inte- 
grated) nominal Ward identity, we arrive at 

d^x J d^y (V;(x){^A^,p}^(x)53(xi)P^(x2)) = 



^ ^ dfs' f d^xi {p\xi)P\x2)) + -4^df^' f d'xi (55(xi)5'(x2)) (70) 
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where the d's are defined in eq. (^). The WI of eq. (^) is vahd if f g and f h, 
otherwise there are additional terms containing flavour-singlet currents on the right 
hand side. Analogous arguments to the ones used previously, show that the ratio 
Zs/Zp, being obtained from a WI, is a finite function of gQ. Thus, the densities S and 
P have the same anomalous dimension. 

Another way of obtaining the ratio ZpjZg is found by combining eqs. ([lOD, (p]), 
(ID and (H): 



Zp_ 

Zs 



niQ — m(mo) 
772-0 - mc 



1 



dWi{g^,mo) 



drriQ 



+ ... 



(71) 



mo=mc 



where uiq — m = 2/9 Z^. 

In this section we have reviewed the lattice vector and axial WI's, emphasizing that 
they explicitly determine finite RC's such as Zy, Za and the ratio Zs/Zp. We have 
repeatedly stressed that these results are strictly true up to vanishing 0{a) terms. In 
the next section we will use the above WI's in order to determine non-perturbatively 
the RC's. A crucial observation is that these computations are performed at fixed finite 
cutoff a. In this case, the continuum limit is never taken, and an apparent "dependence" 
of Zq on a is observed. This "dependence" is a discretization (systematic) error, which 
will eventually drop out in the continuum limit. Minimizing this error is the focal point 



of the Clover action improvement, discussed in subsection 5.2. 



5 Evaluation of RC's from WI's 

The most accurate method for a non-perturbative determination of RC's is based on 
chiral WI's. Only finite RC's can be extracted from them, such as Zy, Za and the 
ratio Zs/Zp. We will use the general results of the previous section in order to show 
how they can be computed non-perturbatively. 

These RC's are determined by requiring that bare correlation functions, calculated 
numerically on the lattice at finite UV cutoff a~^, when renormalized, obey the WI's 
of the continuum theory. This requirement, which we have justified in the previous 
section in the a — > limit, should be also valid at very small but finite lattice spacing, 
within a given accuracy. At fixed lattice spacing, correlation functions can be computed 
numerically, providing us with an intrinsically non-perturbative estimate of each term 
of a WI. These WI's can then be solved for the several RC's. Clearly, this method is free 
from the errors which affect the perturbative estimates of the RC's (i.e. higher orders 
in go). It suffers however, from finite cutoff effects. We shall address this question in 
detail in the rest of this paper. 



5.1 Wilson WI Estimates of Renormalization Constants 

First we consider the WI computation of the RC of the vector current, from suitable 
ratios of correlation functions. We shall use the following symmetrized form of the 
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conserved current: 

Vf = \[V^{x) + V^{x-^i)] (72) 

Being just a symmetrization of V^, it has the same RC, Zyc = 1. From the discussion 
of the vector WI's between hadronic correlation functions, we conclude that Zy can 
be derived from ratios of n-point correlation functions; e.g. 

^ ^ JdxJdye^'^^{PKmo^ix)Pky)) ^ (K(g)|l^o^(0)|Z)(0)) ^ 
' JdxJdye^'^{PKmVo{x)PUy)) " {K{q)\VomD{6)) 



!dx!dye^''-{PKm^{x)Pl{y)) ^ {K{q)\V^mD{^)) 
JdxJdye^'^-{PKmk{x)Pliy)) ^ (K(^)|Ffc(0)|Z)(0)) 



^ ^ JdxJdye^'^-{M>^,{0)V,^\x)PUy)) , Er 4mK*{q)\Vk'' mPjO)) , 
/di;/dyV?-(M^..(0)yfc(x)Pi,(y)) " e^(9l(i^;(^)|14(0)|I?(0)) 

where A; = 1,2,3 is a spatial index and e^{q) (A = 0, ... ,3) the vector meson polar- 
ization vectors. For definitiveness, we have specified the quark flavours by consider- 
ing the correlations characteristic of the p decay (2-point functions) and those of the 
D ^ K and D K* decays (3-point functions). Thus, the operators are defined 
as: Pk = u'-f^s; Pd = u'y^c; = u'-f^d; M^* = uj^s. The vector currents have a 
u — d flavour structure in Rp and a u — c flavour structure in all other ratios. In the 
above equation, we also show the asymptotic behaviour of these ratios, obtained at 
large time separations, up to contributions from higher excited states. All these ratios 
give estimates of Zy up to terms of 0{a). 

The above ratios represent the most straightforward way to obtain Zy, but it is 
by no means the only one. As an example, we present a method, which follows more 
closely the vector WFs but is analogous to the determination of Zy from the ratio Rq. 
It consists in integrating over all 3-space x both sides of eq. (50). Upon integrating, 
the spatial derivative V^V^ vanishes (either because fields die-off at infinity or due to 
periodic boundary conditions of a finite lattice). Then we write the WI in terms of 

rC 




Vn^. The result IS 



dx{P'\xi)V%^{x)P''\x2)) = ^{m2-m3) J dx{P'\xi)[S^Hx)+S^^^{x-b)]P^Hx2)) 

(74) 

where, upon passing from Vq to Vg^, we end up with the sum of two adjacent scalar 
densities on the r.h.s. Finally, we substitute ZyVo for V/f in the above and solve for 
Zy. The implementation of eqs. (73) for the computation of Zy will be referred to as 
the "Ratio" determination. When WI's are directly used instead, we will call it the 
WI determination. Both are non-perturbative methods, equivalent in the a — > limit. 

For the WI determination of Za we start from eq. ([69|). For definitiveness, we will 
consider 3 quark flavours, which may or may not be degenerate in mass. It is necessary 
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to keep track of the flavour content of each operator, e.g. OP(x) = ^i(x)r^/;2(a;). The 
RC's do not depend on flavour. We shah however leave flavour indices on them in 
order to remind us the origin of possible discretization errors in their determination 
(e.g. 0{ami, 01712) for Z^^). With this flavour content, setting X2 = 0, eq. (p9|) becomes 



2p'^j dx j dxi{P^\x)Af{xi)V^\Q)) = 

r 7I3 /• 

= / dMV.''{xi)V^'H0)) - -^r^ / dxi{Af{x,)Afm (75) 

This WI can be useful in several ways. To simplify matters, in practice we will consider 
the case of flavours 2 and 3 degenerate in mass. This means that the pre-factors of 
the two integrals of the r.h.s. of eq. (75) simplify to Z'y / (Z^)"^ and XjZ'y . Thus, 
by solving the above WI at different times ti, we obtain estimates of Z'y and Z^^ . 
Alternatively, Zy can be calculated from, say, one of the ratios of eq. (^) and used in 
eq. (^), from which Za can be computed. Moreover, with 1/ = /? = 0, the integrated 
correlation / (ixi(Vo(xi)Vb(O)) vanishes, since it is proportional to the charge of the 
vacuum. This means that the remaining WI 

2pi2| I dx^{P^\x)Al{x{)Vi\^)) = -^J dx,{Af{x^)Al\0)) (76) 

can be used for yet another determination of Zy- Besides Za_ and Zy, the finite ratio 
Zs/Zp can also be obtained from the WI of eq. ([70|), by a method analogous to the 



one used for Za', eq. (75). 

So far only WI's of hadronic correlation functions have been examined. Quark 
correlation functions such as {ijj{x)V^{0)'ip{y)) and {ip{x)A^(0)il>{y)) can also be used 
in principle with WI's such as eqs. ( p2|) and (|57|). However, in numerical simulations 
the ratios of hadronic correlation functions turn out to be more stable than those of 
quark correlation functions. For the latter case, we must work in a fixed gauge and the 
Gribov ambiguity causes increased statistical fluctuations. Thus, the determination of 
Zy and Za from quark state correlation functions is expected to be noisier. This has 



been explicitly verified for Za in refs. [17| and 



5.2 Clover-Improved WI Estimates of Renormalization 
Constants 

As already pointed out, the WI determination of the finite RC's has the advantage 
of being non-perturbative. In principle, all WI's should give the same finite Zq of a 
given operator O, since the RC's are characteristic of the operator they renormalize, 
but are independent of the ME's from which they are obtained. In practice, however, 
problems arise from the fact that the RC's, properly defined at vanishing cutoff, when 
calculated from WI's (or Ratios), are obtained at finite cutoff. At finite lattice spacing, 
operators mix with higher dimensional (irrelevant) operators. This mixing spoils the 
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renormalization properties of the lattice operators, which is only recovered in the limit 
a — > 0. The contamination from these higher dimensional operators to the numerical 
calculation of RC's depends on the correlation functions (or ME's) used. Consequently, 
the result for Zq obtained at finite cutoff will depend on the correlation functions from 
which it has been obtained and the quark mass at which the simulation was performed. 
This unwanted dependence signals the presence of systematic errors which are effec- 
tively 0(a) in the scaling limit. More specifically, they are ©(oAqc^)), ©(am), 0(ag) 
and higher, where m and q stand for the masses and spatial momenta that characterize 
the process under study. 

In order to reduce this systematic error, several proposals have been put forward 
[51 1, 1 30 1, (3^, [^2|| . All are based on the Symanzik improvement programme [52|, 



which consists in adding suitable irrelevant dimension-5 operators to the lattice action. 
Moreover, d-dimensional operators must also be redefined by the addition of {d + 1)- 
dimensional improving operators. All proposals cited above achieve the elimination of 
the 0{a) discretization error (although at different orders in PT), and are referred to 
as C(a) improvement. The most widely used improved action is the SW/Clover one 
30i, defined as follows: 



a 



X,flU 



lXvF^,y{x)lp{x) 



(77) 



with F^u{x) the clover-leaf discretization of the field tensor. The improving coefficient 



cswido) can be calculated in PT |53]; more recently a non-perturbative determination 
has also been performed |^2|. In this work we limit ourselves to its tree- level value 
csw = 1- At this order, all 0{agQ^ln" a) terms, which are effectively of 0{a) in the 
scaling limit (qq ~ 1/lna), are eliminated from correlation functions. At leading- log 
level, the improvement of local operators can be expressed as a rotation of the fermion 
fields pT|. According to this prescription, improved bilinear operators are given by 



Of (x) = ij^{x)Tij''{x) 



(78) 



where the rotated fields are defined through 



tP^ix) 
'ip^{x) = -ipix) 



!--(]/) -mo 



ipix) 



1 + - ( ^ +mo 



(79) 



and the symmetric lattice covariant derivatives are: 

1 



a -tpix) = -[U^{x)^{x + A) - Ul{x - fi)ip{x - fi)] 



+ f^)Ul{x) - i:{x - fi)U^{x - fi)] 



(80) 



The improved operators differ from the original ones by terms proportional to the 
cutoff. Consequently, they have different RC's, say Zqi. 
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For the Clover action, the quark propagator S^{x — y) = {■ip{x)ip{y)) satisfies the 
equations 



p {x) + mo + aW 



S\x -y) = 5{x - y) 



S\x-y) 



p {y) + mQ + aW 



5{x - y) 



(81) 



The terms W and W of eq. ( |8l| ) are shorthand notation for the operators arising from 
the variation with respect to ip and ip of the Wilson term of the fermionic action, 
including the Clover term given by eq. (77). The tree-level 0{a) improved quark 
propagator is {ip^{x)'4>^{y)) . Using eqs. (79) and (|8l[), we can write it in terms of 
S'{x-y): 



(V^^(x)V^^(y)) 



Ri 



1 



y) 



+ -6{x 



y) + 0{a^) (82) 



where 0{a^) terms proportional to W and W have been dropped, as they do not 
affect 0{a) improvement. In computing correlation functions between external on-shell 
hadron states, the ^-function never contributes. Therefore these correlation functions 
can be expressed directly in terms of the effective rotated propagator 



S'ff{x-y) 



1 



S\x - y) 



(83) 



This is the Clover propagator we will mostly use in our computations. The reason it 
contains an 0{a?') term in its definition is that it is readily computable in numerical 
simulations |T^. Thus, it has been extensively used in several Clover improved lattice 
QCD computations of on-shell ME's. However, upon studying the lattice WI's we find 
cases for which the space-time points of correlation functions are allowed to coincide. 
The on-shell ME's argument is then invalidated and the complete propagator of eq. ( ^2|) 
must be used, with the (5-function giving rise to contact terms. 

More generally, we can add to the improved propagator ('0^(x)V^^(y)) the two 
equations of motion (^TJ), multiplied by an arbitrary factor z. Up to 0(0^), we then 
find 



(V'«(x)V^^(y)) 



1 



z ]p [x) — (1 — z)mQ 



S\x-y) 



l + -[P{y)z + {l 



z)mo 



+-{2z-mx-y) + 0{a') 



(84) 



In this case, improved on-shell correlation functions can be expressed in terms of traces 
of the effective propagator 



SfHx-y) 



1 



-(zp{x)-{l- z)mo 



S'{x-y) 



1 + - ( ^ {y)z + (1 - z)mo 



(85) 
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with z a free parameter, which we have at our disposal for optimization purposes. 
Clearly, for different choices of z, these ME's (and the corresponding improved oper- 
ators) differ by ©(a^) terms. This implies that their RC's also vary with z. In this 
work, unless otherwise stated, we shall use operators Of which correspond to z = 1. 

Although on-shell correlation functions of local improved operators Of are now im- 
proved up to C'(g'oa), the conserved (point-split) current is not improved. An improved 
conserved current is given by 

y^\^) = + \a^p{^a,^i,) (86) 
where the extra "(jp^ term" has a symmetric lattice derivative (2o^ ^f(x) = f{x + fi) — 

Having improved both local and conserved operators, we can now use them in the 
ratios of eq. (^), in order to obtain estimates of Zyi which only suffer from ©(a^g) 
discretization errors. On the other hand, the l.h.s of the WI of eq. (^), even when 
expressed in terms of and Vq"^ , still suffers from 0{a) corrections. This is due to 
the presence of the asymmetric lattice derivative = dl^ + 0{a). In order to improve 
it, we must express it in terms of the symmetric derivative = d!^ + 0{a^) and make 
use of the fact that the conserved currents and V^^ have the same divergence 
(V^V^ = V^V^^). Thus we arrive at the following 0{a) improved version of eq. (^4[): 

dx{P'\xi)vlv,^\x)P^\x2)) = (87) 

_ {m2-ms) f ,-j,„i2/ ^nc23^^^ , c23/^ , n\ , c23/ 



J dx{P^\xi)[2S^^{x) + S^^{x + 6) + S^^ix - 0)]P^\x2)) 



In the above equation we have dropped the I superscript from the density P. The 
factor 1/4 on the r.h.s. comes from V^. We can now express Vg^^ as ZyiV^ and solve 
for Zyi. This estimate of the RC is free of leading 0{a) errors. 

We now turn to the WI of eq. ( |75[) which, expressed in terms of improved operators, 
can be used for the determination of the tree- level 0{a) improved RC Z^^/ of the axial 
current. The subtlety is that the integration / d'^x of the WI requires inclusion of the 
contact terms arising from the (5-function of the rotated propagator of eq. (|84[). Thus, 
for any value of the z-parameter, the Clover improved version of eq. ([75|) becomes 



2p^'j dx J dy{P'\x)AlHy)V^'HO)) = - p'^az) J dy{V^\y)V^\Q)) 

-(wIsT - P''^^) I dy{Af{y)Afm (88) 

where z = {2z — 1). All operators in the above equation, as well as the 2p factors, 
are to be understood as improved quantities. For 2p this is true provided that it is 
computed from eq. (U) with improved operators Aj^ and P^ and a symmetrized lattice 

temporal derivative V*^. 
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5.3 Comparison of the Perturbative and WI Estimates of 
the RenormaUzation Constants 

From the discussion so far, it should be clear that it is essential to the reliability of 
any WI computation, performed at finite cutoff, that any scale with mass dimensions 
(such as Aqcd, masses m and momenta magnitudes q) must satisfy the condition 

AQCD,m,q <.a-^ (89) 

In an ideal situation (g^ and a extremely small) the above condition is satisfied and the 
estimates of the Z's obtained from different correlation functions with the WI/Ratio 
method should all agree. They should also agree with the Z estimates obtained in PT, 
carried out at sufficiently high order. The argument can now be turned upon its head: 
if the Zq^s obtained from WI's or ratios of different correlation functions are not in 
good agreement, this is to be interpreted as a signal that the calculation suffers from 
large discretization errors. In fact, the "dependence" of Zq from the external states 
and/or the quark mass is a criterion used to estimate the finite lattice spacing effects 
which affect the lattice calculations. On the other hand, if the WI/Ratio result for Zq 
turns out to be reliable, any discrepancies with its PT estimate is to be attributed to 
the truncation of the perturbative series. 

To illustrate the point, we collect some results in table ^, obtained at a typical value 
of the lattice coupling constant, (3 = 6.0. For the Wilson case, the presence of large 
0{a) effects is highlighted by the 20% — 30% variation in the various non-perturbative 
estimates of Zy- This wide spread of values renders problematic their comparison to 
the perturbative estimates (both SPT and BPT). For Za, the non-perturbative values 
are compatible within statistical errors. However, these errors are too large to allow 
any conclusion on the accuracy of the two perturbative determinations. 

For Clover fermions, the spread of the different WI estimates is reduced to less 
than 5%. This is due to the elimination of 0{a) terms guaranteed by the Clover 
improvement. This result implies that lattice QCD calculations, performed with the 
Clover action at light quark mass, are significantly less affected by the finiteness of 
the lattice spacing than their Wilson counterparts. Note that BPT predictions for Zy 
agree well with the non-perturbative results. However, for the ratio Zs/Zp BPT is 
less successful. Moreover, the PT estimate of Z^ is less than 1, whereas the NP value 
exceeds unity. BPT cannot rectify this discrepancy. Thus, 0{g^) terms are probably 
important in these cases. 

So far we have presented results at /? = 6.0. The reason for this choice is that this 
coupling is typically used by present-day lattice QCD groups for calibration purposes. 
Most studies of an explorative nature (e.g. RC's from WI's, or from the NP method, 
with or without Clover improvement) have first been carried out at /? = 6.0. Thus a 
comparison of the RC's obtained from all different methods is possible at this value of 



the coupling. In subsection |7.1| we also present results for Zy obtained from the ratio 
method with the Wilson action at /? = 6.4. Although the cutoff is considerably 
increased, we will see that finite cutoff effects remain sizeable. Thus, our qualitative 



27 





Method 


Wilson 


Clover 




SPT 

MC-PB(lnn) 
Ratio [Rp) 
Ratio [Rq) 
WI (Axial) 
NP 


0.83 
0.63 

0.57(1) 1 
0.73(3) 
0.79(4) 
0.75(5) 


11 
l| 

H 


0.90 
0.80 

0.82 [This work] 
0.82 [This work] 
0.80(2) [This work] 
0.84(1) |2§ 


Za 


SPT 

MC-PB(lnn) 
WI (Axial) 
NP 


0.87 
0.72 

0.85(7) 1 
0.80(5) 


B 

241 


0.98 
0.95 

1.11(2) [This work] 
1.06(8) P 


Zs/Zp 


SPT 

MC-PB(lnn) 
WI (Axial) 
NP 


1.10 
1.34 

N/A 
N/A 


1.20 
1.54 

1.64(5) [This work] 
1.6(1) ^ 



Table 3: Estimates of the EC's of the vector current (Zy), the axial current (Za) and the 
ratio of the scalar and pseudoscalar densities (Zs/Zp) from various perturbative and non- 
perturbative methods (see text for labels) at {3 = 6.0. All results are obtained at a small value 
of the light quark masses (am ~ 0.07^. The corresponding references are shown in square 
brackets. 



conclusions on systematic errors characterizing the Wilson action remain valid at /? = 
6.4. Finally, finite RC's with the Clover action (at light quark mass) have also been 



computed at /3 = 6.2 [2C]. 



6 Lattice - Continuum Normalization of Fields 
and Currents 



We now examine a proposal due to Lepage, Mackenzie and Kronfeld [36| -[38] for the 
removal of 0{am) systematic errors which are relevant in lattice calculations with the 
Wilson action at large quark masses. After a critical review of their recipe, we extend 
it to the Clover case. For both the Wilson and Clover cases, we apply it to the WI 
computation of Zy and Za and discuss its merits and shortcomings. 

First we set up our notation. So far we have worked with the lattice action S = 
Sf + Sg of eqs. (|8|) and (^), which depends on lattice fields and bare parameters; i.e. it 
has the form S['ip,ip,Up_,mo, g^, a]. In numerical simulations, it is convenient to work 
with dimensionless fermion fields and the hopping parameter instead of the quark mass; 
i.e. the action used in Monte Carlo simulations 



gLAT 
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+ T.x{x)x{x) + Sg (90) 

X 

has the form S^^^[x, Xj U^, K, Qq], where K and x ^-re defined through 

X = ]l^^ , ^-4 = a-o (91) 

Note that the lattice spacing a has been scaled away in S^^"^. All the equations 
appearing in the previous sections can now be trivially rewritten in terms of x a-^d X 
fields and masses in lattice units moa (or K); the dependence on a drops out and the 
equations express relations between dimensionless quantities. From now on, all lattice 
quantities (such as quark propagators, operators and their correlation functions) are to 
be understood as defined in terms of the x and x fields; i.e. they are (dimensionless) 
quantities expressed in lattice units. 

In order to connect lattice fields to continuum (bare) ones, we must expand the 
lattice action for a ^ at fixed go and K. This action then reduces to the continuum 
QCD action, S^^^ , provided we relate the continuum field Af^ to its lattice counterpart 
Ufj^ in the following way: 

Uf, = exp{iagoA^) (92) 
The standard prescription for the transcription of any continuum quantity on the 



lattice follows from eqs. ( |9l|) and (92). For example, in the naive continuum limit, 
the dimensionless lattice quark propagator S{x — y) = {xix)x{x)) becomes its QCD 
counterpart S^^^{x — y) = {ip{x)il){x)) when normalized as follows: 

a'^ 2K S{x-y) ^ S^^^{x-y) (93) 

Clearly, such transcriptions are not unique. Modifications of the lattice quantity by 

terms which are 0{a) give equally acceptable transcriptions. For example, both the 

lattice local vector current V^(x) = Xi{x)l ij.X2{x) and the corresponding point-split 

conserved current {x) are good transcriptions of the vector current obtained from 
gQCD._ 



a 



-3 



^2Ki2K2V^{x) ^ K^^^(x) 



a-V2i^i2E:2 y^^(x) ^ y^^^^(x) (94) 

Note that these relations are just transcriptions of bare QCD quantities in terms of 
quantities derived from S^^^ . They are not relationships between lattice and renor- 
malized quantities, which must involve the presence of RC's; e.g. Vf^ = ZyV^. In our 
notation, the QCD superscript and the "hat" distinguish the two cases. 



6.1 KLM Factors for the Wilson Action 

The proposal by Lepage, Mackenzie and Kronfeld (KLM) [^6| -||3^ is a refinement 
of the above matching. Their aim is to account for large 0{am) corrections in the 
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limit of small aq {q stands for the magnitude of any relevant spatial momentum). 
In practice their recipe consists in calculating the matching factor between lattice 
correlation functions and their continuum counterpart in the free theory, for finite 
lattice spacing and vanishing spatial momentum. The claim is that this free theory 
correction should account for most of the 0{am) systematic error also in the interacting 
case. For example, in the free theory, the matching between the continuum and lattice 
quark propagators in (t,p)-space at p = is 

S^'^^{t,p = 0) = a-^2K{l + amo)S{t,p = 0) (95) 

This implies for the lattice fermion field the matching relation 



^1/2 [2K 



2K , , 

3 X (96) 



a 

where the KLM wave function normalization factor is 

F^iam^) = 1 + amo (97) 

This factor tends to 1 in the limit amg — > and the standard matching of eq. (^) is 
recovered. We stress that this is a free theory matching which is due to the granularity 
of the space-time. It is not a renormalization. 

A similar procedure applies to all correlation functions. For example consider the 
correlation function of the lattice local current between external quark states at zero 
spatial momentum, 

L^it^.ty) = E(xi(y)^/.(^)X2(0)) (98) 

It can easily be seen that, in the free theory and at finite cutoff, it is related to the 
continuum correlation function by 

L^^"{tx,ty) = a-^2Ki2K2F^{amoi)F^{amo2)Lf,{U,ty) (99) 

where aniQi and amo2 are the bare quark masses of the two distinct flavours of the 
quark fields. Thus, there is no extra KLM normalization, apart from the one related 
to the quark fields. The KLM factor for the local vector current turns out to be 1. 
This is also the case of the local axial current and of any other local bilinear operator. 
The KLM factor for the spatial component of the conserved current V^'{x) is also 1. 



On the other hand, as pointed out in |54], the temporal component of the conserved 
current Vq'{x) requires an extra KLM factor 

Fyc{amoi,a'mo2) = ( H ^ ) (100) 

It arises from the matching of 

Co{t,,ty) = J2ixiiy)Vo''{x)x2m (101) 
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to Cq^^ {ty,tx)i which gives: 

Cq {ty,tx) = a~ 2Ki2K2F^{amoi)F^{amo2)Fyc{amoi,amo2)Co{ty,t^) (102) 

It is now trivial to derive the KLM matching factors of hadronic correlation func- 
tions. We just have to consider the number of quark fields and the composite operators 
entering in the correlation function. For example, for the 2-point function of the vector 
currents, with the /?-meson quantum numbers, we have 

X X 

whereas for the vector 3-point function of the D ^ K decay we have 

^(PK(0)yo^(x)pt,(y))«^^ = (104) 



= a''^2Ku2Ks2KcF^{amou)F^{amos)F^{arnoc)Fyciarnos,arnoc) ■ 

The KLM factors obtained above depend on the bare quark mass of eq. (H). A 
possible improvement, which could take into account interaction effects, consists in 
using a more physical definition of the quark mass in the KLM factors. A natural 
choice is, for instance the subtracted quark mass m which, in terms of the hopping 
parameter K, is given by 

am = (105) 

2K 2Kc 

with the critical value Kq determined non-perturbatively. Unless otherwise stated, all 
KLM factors used in our numerical analysis are calculated at the above value of the 
quark mass (e.g. F^{am), Fyc{a'm) etc.). Another proposal has been made in p9| , 
based on the observation that in lattice perturbation theory, a large fraction of renor- 



malization effects comes from tadpole diagrams. The authors of |2S] propose to take 
the bulk of these effects into account by implementing a Mean Field Tadpole Improved 
(MFTI) prescription. MFT Improvement amounts to the following two substitutions: 



'1) The link variable in any operator, such as Kp, is to be substituted by 



no 

where no is any reasonable mean-field estimate of the expectation value of the link. 
Two standard choices are no = [^TrUpY^^ and no = l/{SKc). In practice, for current 
values of the coupling constant, the two no estimates differ by about 10% (about 2% 
for the Clover action). In the results of this work, the latter estimate of no has been 
used. 

(2) The hopping parameter is obtained by substituting 

k = Kuo (107) 
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This implies, for example, that the MFTI quark mass is obtained by substituting 

amo ^ am = - 4 = ^Kc{^ - ^) = SKcam (108) 
2a ^a 



The "standard" MFTI normalization of the quark field is given by the factor v 2K; its 
KLM/MFTI correction is obtained through the substitution 

V2KE^{am) \[2kF^{am) = \j2k{l + am) (109) 
while for the conserved current Vq , the KLM/MFTI factor becomes 

It must be kept in mind that the 1/tio factor comes from the link C/^ present in the 
conserved current (i.e. it is a redefinition of the lattice operator due to the redefinition 
of the gauge field through eq. ( [106D ) whereas the mass dependent factor arises from 
eq. ( |108D (i.e. it is a rescaling of the bare mass). This difference in the origin of 
the two factors will be important when we compute ME's of the conserved currents 



with KLM/MFTI corrections (see comments after eq. (124)). We also note that the 



prescription of eq. (|108|) has also been obtained by a tadpole resummation argument, 



without recourse to MFTI, in ref. [42|. 



The KLM factors, having been calculated in the free theory and at zero quark spatial 
momentum, do not account for, say, 0{gQam), 0{aKQCD) and 0{ap) corrections. The 
hope is that at large quark mass and small coupling, they account for most of the 
systematic error due to the finiteness of the UV cutoff. In sect. |^ we will test these 
expectations. 



6.2 KLM Factors for the Clover Action 



We now calculate the KLM factors for the Clover case. Since the free fermion prop- 
agator is the same for both the Wilson and Clover actions, its KLM factor is the 
same. However, fermion fields need to be rotated, as for example in eq. (l79|). Since we 
usually opt for this specific rotation (z = 1), most correlation functions are expressed 
in terms of the effective lattice propagator S^^^ [x — y), given by eq. (p3[). The KLM 
factor for the effective propagator is then given by 



5«^^(t,p = 0) = a- 



2K{l + amo) 



1 

1+4 



1 + amo 



1 



■ S^ff{t,p = 6) (111) 



1 + amo , 

where again everything is worked out in the free case. Thus, compared to the Wilson 
case of eq. (p5|), there is an extra term, Fji{amo), arising from the fermion rotations in 
the improved KLM factor: 



Fi(amo) = F^{amo)FR{amo) = (1 + amo) Fji{amo) 



(112) 
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with 



FR{amo) 



1 + 7(1 + amo 



-2 



(113) 



4 V 1 + aiTio , 

It can be seen immediately that the overall KLM factor, F^laniQ), is 1 + 0{amQ), in 
accordance with 0{a) improvement. From now on, we will directly substitute the bare 
mass mo by the subtracted one m, in the spirit of eq. ( |105D . 

Let us consider the quark correlation function of the local improved current = 
Xi7fiX2 (i-®- improved analogue of eq. 



Li{t.,ty) = E(xf (y)^/(^)xl'(o)) (114) 

It is matched to its continuum counterpart through 

^QCD ^ a-'F^{am,)F^{am2)L'^{U,ty) (115) 

Thus we see that also in the Clover case the KLM factor for the local vector current, 
as well as for any other local bilinear operator, is 1. 

We now pass to the KLM factors of the conserved improved current of eq. (p6[). 
Unlike the improved local current V^, fermion fields are not rotated for this current. 
We first consider the improved quark correlation function of the current's temporal 
component Vq^^: 

The matching factors for this correlation are easily found to be: 

C^''''{ty,t,)= (117) 

fi 1/2 1/2 T 

a~°2Ki2K2F^{ami)F^{am2)Fj^ {ami)Fj^ {am2)Fyc{ami , arn,2)Co (ty, tx) 



where Fyc is given in eq. (|10C| ). Note that at p = the Vi_iX^fiOX term of the current 
vanishes. The above overall KLM factor in eq. ( |117| ) is 1 + 0{a^m?)^ in accordance 
with improvement. 

Finally, we examine the improved correlation Cl{tx,ty) of the spatial component 
of the conserved improved current: 

ci{ty,tx) = Y.^xf{y)v^\x)xSm (118) 

x,y 

The standard KLM computation yields the following matching: 

C'^''''{ty,t.)= (119) 

(\ 1 /2 1 /2 7 

a~ 2Ki2K2F^{ami)F^{am2)Fjl {ami)Fj^ {am2)Fyci{ami,am2)Cf^{ty,tx) 
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where 



1 



1 — (1 + ami)^(l + am2) 
4(1 + ami)(l + ani2) 



1 -1 



(120) 



Again, being 1 + 0{am?), the overall KLM factor in eq. ( 119 ) is consistent with im- 
provement. 

We conclude by noting that using the Clover version of the axial WI, in the spirit 
of KLM, requires a matching factor also for 2p. From the improved version of eq. (|64[), 
with a symmetric lattice temporal derivative v'^, we obtain the KLM matching 



2p' 



QCD 



2{ami + am2) 



exp(ami + om2) — exp(— ami — o,m2 



-2/ 



(121) 



Just like in the Wilson case, m's and mq's can be introduced, whenever required 
by the prescription, in order to account for KLM/MFTI matching factors. Again, 
the extension of the Clover KLM factors to those characterizing hadronic correlation 
functions is straightforward. 



7 Non-perturbative Results for the RC's 

In this section we present results for the RC's Zy and Za and the ratio Zs/Zp obtained, 
non-perturbatively, from numerical studies of the chiral WI's on the lattice. The RC 
of the local vector current, Zy, has been also computed by taking the ratio of ME's of 
the conserved to the local vector current between different external states. 

The main purpose of this study is to investigate the systematic errors present in 
these calculations and, in particular, those due to the finiteness of the lattice cutoff. 
For this purpose, we have considered both the standard Wilson action and the tree- 
level C'(a)-improved Clover action. In order to enhance the finite cutoff effects, we 
have studied a large range of values of the lattice bare quark masses, which vary from 
am ~ 0.05 up to am ~ 1. As expected, we find that these effects are typically larger 
with the Wilson action than with the Clover action. 

Moreover, we have also tried to correct the finite cutoff effects by following the 
KLM prescriptions, discussed in the previous section. We find that KLM corrections 
are typically more effective in the Wilson case, where they are 0{am) rather than in the 
Clover case, where they are O^a^m^). The analysis of the Clover RC's shows that the 
bulk of the finite cutoff effects depends linearly on the bare quark mass. This suggests 
the existence of large 0{g^am) terms, which, being a consequence of the quantum 
loops, cannot be corrected by the tree-level KLM prescriptions. 

A significant result is that the finite cutoff effects, even in the Clover case, can be 
as large as 20 — 30%, for values of the bare quark masses which are currently used in 
numerical simulations (e.g. am ~ 0.5). This result strongly supports the necessity of 
an additional systematic improving strategy, such as |^, In the following we 

present our results in detail. 
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7.1 Ratio Determination of Zy with the Wilson Action 

We first present a computation of tlie RC of the local vector current, Zy ^ in the 
case of the standard Wilson action. The data have been obtained by averaging over 20 
quenched gauge field configurations, generated at /3 = 6.4. This value of the lattice cou- 
pling constant corresponds to a relatively large value of the lattice cutoff, a^^ ~ 3.6—3.8 
GeV. The lattice volume is x T = 24^ x 60. The same set of gauge configurations 
and quark propagators has already been used in previous phenomenological studies of 



heavy meson physics on the lattice |55]. 



Zy has been computed by considering different ratios of two- and three-point corre- 
lation functions which reduce, asymptotically in time, to different ME's of the conserved 
and the local vector currents; cf. eq. (|73|). For simplicity, these ME's will be denoted in 
the following as < 0|Vfc|/9 >, < K\y^j\D > and < K*\y^\D >, regardless of the specific 
flavour content of the mesons in the external states. In particular, by D we mean a 
pseudoscalar meson consisting of a "light" and a "heavy" quark, whereas by i^, K* and 
p we indicate mesons composed by two light quarks. In all cases, the Wilson parameter 
of the light quark mass is fixed at K = 0.1485, the critical value being Kq = 0.1506. 
Thus this mass corresponds to a bare quark mass am ~ 0.05. The heavy quark mass 
varies in the range between am ~ 0.2 and am ~ 0.6. In the three point correlation 
functions, the source of the final K and K* mesons is inserted at the time t = 0. The 
D meson is at t = 28 and the vector current varies in the interval 12 < t < 16. The 
D meson is always considered in its rest frame and the spatial momentum carried by 
the current is denoted by q. For simplicity, the momentum transfers q = (0, 0, 0) and 
q = (1,0,0) (expressed in lattice units of 2tt / La) will be abbreviated as g = and 
q = 1 respectively. 

A first way of obtaining Zy consists in computing the ratio Rq of three point 
correlation functions defined in eq. ([73|). This ratio has been computed at momentum 
transfer q = In the large time limit, it is given by: 

According to the KLM claim, the bulk of the 0{a) corrections in this expression should 
be given by the KLM factor Fyc{ami,am2) of eq. ( |100D associated with the vertex 

of the conserved current Vq" (the KLM factors F^{am), associated with the external 
quark fields, cancel in Rq). Therefore we also define the ratio: 

= R, . Fyc{ami,am2) = Ro ■ U + ^!!!l±^] = Zy + 0{a) (123) 



If the KLM prescription works in this case, the 0{a) terms in eq. (123) should be 
smaller than those in eq. ( |122D . This means that the ratio Rq^^^ should exhibit a 
milder dependence on the two quark masses. In addition, we also consider the MFTI 



version of the ratio Rq. It is obtained from eq. ( |123|) , by substituting the bare quark 
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K2 


anieff 


Ro 


nKLM 

Kq 


Kq 


0.1275 


0.1275 


0.602 


1.278 


0.798 


0.741 


0.1275 


0.1325 


0.528 


1.207 


0.790 


0.738 


0.1275 


0.1375 


0.459 


1.141 


0.782 


0.735 


0.1325 


0.1325 


0.454 


1.136 


0.782 


0.735 


0.1275 


0.1425 


0.395 


1.081 


0.775 


0.732 


0.1325 


0.1375 


0.385 


1.070 


0.773 


0.731 


0.1275 


0.1485 


0.324 


1.017 


0.768 


0.731 


0.1325 


0.1425 


0.321 


1.010 


0.764 


0.728 


0.1375 


0.1375 


0.316 


1.005 


0.763 


0.728 


0.1375 


0.1425 


0.253 


0.945 


0.754 


0.725 


0.1325 


0.1485 


0.250 


0.946 


0.756 


0.727 


0.1425 


0.1425 


0.189 


0.885 


0.744 


0.721 


0.1375 


0.1485 


0.182 


0.880 


0.745 


0.722 


0.1425 


0.1485 


0.118 


0.820 


0.733 


0.718 



Table 4: The ratios Rq, Rq^^^ and Rq^^'^^ for the Wilson action and for different values of 
the quark masses. In the three point correlation functions the light spectator quark has mass 
0.05. The statistical errors in the table are less than (or equal to) 1 on the last digit. 



am 



masses with their MFTI expressions: 

/ ami + arh2\~'^ 
^ 2 ) 



T>MFTI _ jj 



Zv + 0{a) 



(124) 



Notice that, with respect to the definition of the KLM/MFTI factor Fyc given in 



eq. (110), we have omitted in the previous equation a factor uo- As pointed out after 
eq. (lie), its presence would imply a redefinition of the conserved current, whereas 
we are only interested in a comparison of the discretization errors with and without 
KLM/MFTI corrections. Therefore, by omitting uq, we are computing the ratio Rq^^'^^ 
of the ME's of the same operators which appear in the two ratios Ro and Rq^^^. In 
this way, all three ratios give estimates of the same Zy and only differ by 0{a) effects. 

The values of the three ratios i?o are shown in table ^ for different values of the 
hopping parameter. They are also plotted in figure |l| as a function of the "effective" 
quark mass, am^ff = {ami + am2)/2. We see that Rq has a strong dependence on the 
bare quark masses, varying by approximately 50% as the mass varies from am ~ 0.2 to 
am ~ 0.6. This dependence is however significantly reduced when we apply the KLM 
correction. This reduction is somewhat bigger in the MFTI prescription. 

In order to discuss these results in a more quantitative way, we fit the three ratios 
Rq, asa. function of the quark masses, taking into account contributions up to quadratic 
terms: 



Rq = A + B{ami + 07712) + C{ami + am2f' + D{ami — am-2)^ 



(125) 
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Figure 1: The various ratios Rq, defined in eqs. i\12^ ), (\12^ and ( \124\ ) for Wilson fermions, 
as functions of the effective quark mass. 



Since the current vertex is symmetric under interchange of the two quark fields, there 
is no linear dependence on the mass difference (ami — ara2). For all three ratios i?0! 
Rq^^ and Rq^^^^ the value of the coefficient A in eq. (|125| ) gives an estimate of Zy, 
with presumably negligible 0{am) systematic errors, which we denote by Zy^. The 
result in all three cases is approximately the same and we quote, as an overall estimate: 

Zl/^ = 0.71(1) , from i?o (126) 



The parameter B in eq. ( |125D gives the size of the leading 0{am) corrections. We find 
that B decreases from B ~ 0.46 for Rq, to ~ 0.09 for Rq^^^ and to an even smaller 
value, B ~ 0.03, for R^^^'^^ . The coefficients C and D of the quadratic terms are 
always of the order 10~^ or smaller in magnitude. Thus, in the case of the ratio Rq, 
the KLM prescription correctly takes into account and removes the bulk of the finite 
cutoff effects. 

Next we compare results for Zy obtained from different ratios of correlation func- 
tions, cf. eq. ([73[). In particular, for different values of the quark masses, we have 
computed: 



^p{Q = 0) = — , — ^ - Zy + U[a} 

Erer(0|Vfc(0)|p/ 



^'^"^ - ^) - {K\vm\D) - + ^^"^ ^^^^^ 
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Ki 


K2 


arrieff 


Rp 


Rk 


^k 


0.1275 


0.1275 


0.602 


0.622(1) 


1.073(41) 


0.643(23) 


0.1275 


0.1325 


0.528 


0.628(1) 


0.972(37) 


0.639(23) 


0.1275 


0.1375 


0.459 


0.633(1) 


0.877(34) 


0.634(23) 


0.1325 


0.1325 


0.454 


0.632(1) 


0.979(38) 


0.634(23) 


0.1275 


0.1425 


0.395 


0.638(1) 


0.792(30) 


0.628(22) 


0.1325 


0.1375 


0.385 


0.637(1) 


0.883(34) 


0.630(23) 


0.1275 


0.1485 


0.324 


0.643(2) 


0.719(20) 


0.625(21) 


0.1325 


0.1425 


0.321 


0.642(1) 


0.795(30) 


0.624(23) 


0.1375 


0.1375 


0.316 


0.641(1) 


0.891(34) 


0.625(23) 


0.1375 


0.1425 


0.253 


0.645(1) 


0.801(30) 


0.620(23) 


0.1325 


0.1485 


0.250 


0.646(2) 


0.721(20) 


0.623(22) 


0.1425 


0.1425 


0.189 


0.649(1) 


0.812(30) 


0.616(23) 


0.1375 


0.1485 


0.182 


0.649(2) 


0.723(19) 


0.621(22) 


0.1425 


0.1485 


0.118 


0.653(2) 


0.727(19) 


0.618(23) 


0.1485 


0.1485 


0.047 


0.658(2) 







Table 5: The ratios Rp, Rk and Rl for the Wilson action at different values of the vector 
current quark masses. In the three point correlation functions the light spectator quark has 
mass am ~ 0.05. 



The KLM corrections to the correlation functions of these ratios always cancel out 
between numerator and denominator. Thus, according to the KLM expectations, all 
these ratios should be free of large 0{am) corrections. 

The values of Rp, R^ and R^. are given in table ^ for different values of the vector 
current quark masses. The same ratios are also plotted in figure ^ as a function of the 
effective quark mass. In the figure, the ratio R^^^^ is also plotted for comparison. 
We see that Rp, Rk and R^. exhibit quite a different behaviour as a function of the 
quark masses. When the mass varies from am ~ 0.2 to am ~ 0.6, the ratios Rp 
and R^ remain almost constant, while the ratio Rt is found to be affected by large 
0{a) corrections. In particular, by performing a fit as in eq. ( |125[ ), we find a large 
contribution coming from the term proportional to (ami — am2)^, with a coefficient 
D ~ —0.51. Compared to its values corresponding to the three previous estimates of 
-Ro, D of Rk is larger by approximately an order of magnitude. The total finite cutoff 
effects in the ratio Rk can be as large as 50% for am ~ 0.6. These effects cannot be 
corrected by the KLM prescription. 

The values of the coefficients A of the fits give an estimate of Zy with presumably 
negligible 0{am) effects. We quote: 





= 0.66(1) 


from 


Rp 






= 0.68(1) 


from 


Rk 


(128) 


Zy 


= 0.62(1) 


from 


^k 
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ure 2: The ratios Rp, R^^^^ , R^ and RI, for Wilson fermions, as functions of the effective 
\rk mass. 



The differences between tliese values and tfie Rq estimate of eq. ( |126[ ) for Zy represent 
a measure of the 0{a) corrections which remain after the 0{am) terms have been 
presumably removed by the fit. These corrections, which are expected to be 0{aq) or 
0{aKQCD)i leave a systematic uncertainty which is approximately 15%. 

7.2 Ratio Determination of Zy, Za and Zp/Zs with the 
Clover Action and Light Quark Masses 

In this subsection we present the results of a calculation of the RC's Zyi ^ Zj^i and the 
ratio Zpi /Zgi obtained with the tree-level C'(a)-improved Clover action, by using the 
lattice chiral WI's. From now on, the superscript "/" of the Zqi^s will be dropped. We 
limit ourselves to light values of the bare quark masses. The value of the lattice coupling 
constant is /? = 6.0. The lattice volume is 16^ x 32. We have generated 18 quenched 
gauge field configurations. We have worked at three values of the Wilson hopping 
parameters, namely K = 0.1410, 0.1425 and 0.1432 corresponding to bare quark masses 
ranging from am = 0.056 to am = 0.111. The critical value is Kq = 0.14554. 

All two- and three-point correlation functions have been computed by considering 
degenerate quark masses. The results are a replication of those of ref. ||l^ for three 
quark masses. The RC of the local axial current Za and the ratio Zp jZg have been 
obtained from gauge invariant chiral WI's, between external hadron states (see eqs. ( |69| ) 
and (|70|)). Gauge dependent chiral WI's between external quark states have also been 
implemented for a calculation of Zj^ (see eq. (^)). The RC of the local vector current, 
Zy ^ has been computed from the ratio i?p, with zero spatial momenta. The results of 
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Ki 


ami 


Zv 


Za (g.i.) 


Za (g.d.) 


Zp/Zs 


0.1410 
0.1425 
0.1432 


0.111 
0.073 
0.056 


0.833(1) 
0.824(2) 
0.819(2) 


1.06(5) 
1.09(3) 
1.07(7) 


1.07(7) 
1.14(8) 
1.17(13) 


0.61(3) 
0.60(2) 
0.60(3) 



Table 6: The RC's Zy, Za and the ratio Zs/Zp for the Clover action with light quark 
masses. The labels "g.i." and "g.d" distinguish results obtained from gauge-invariant and 
gauge-dependent WI's respectively. 



this computation are presented in table g. 

We see that, in the case of Zy, the statistical errors are sufficiently small to reveal 
the presence of finite cutoff effects in the calculation. By performing a fit as in eq. (125) 
we can obtain the value of Zy, with presumably negligible am effects. We quote: 

Z^^f = 0.80(1) , from Rp (129) 

As seen in table the variations of the RC Za and the ratio Zs/Zp with the light 
quark masses are of the size of the statistical errors. Thus, in these cases the presence 
of finite cutoff effects cannot be determined. An analogous study at /3 = 6.2 has been 
performed by the UKQCD collaboration |2C]. 



7.3 Ratio Determination of Zy with the Clover Action 
and Heavy Quark Masses 

In order to investigate the size of finite cutoff effects when the tree- level C'(a)-improved 
Clover action is used, we have also computed Zy at several large values of the quark 
mass. We have generated 30 quenched gauge field configurations at /3 = 6.0. The 
lattice volume is 16^ x 48. Two different values of the light quark hopping parameter 
have been considered {K = 0.1425 and K = 0.1432) and four different values of the 
heavy quark mass. The latter correspond to bare quark masses varying from am ~ 0.3 
up to am ~ 0.9. 

We have computed Zy from the ratios i?o and Rp of eq. (|73|), with quark propagators 
improved according to eq. (|83|). All external states and currents are always at zero 
spatial momentum. In the three-point correlation functions the pseudoscalar D meson 
is inserted in the middle of the lattice, at i = 24, so that results from both positive 
and negative times can be averaged, thus improving the statistics. 

In order to investigate the effects of the KLM corrections in the Clover case, we 
have also considered the corresponding KLM-improved versions of the ratios Rq and 



Rp. According to eqs. ( |115|) , (|117| ) and (|119| ), these corrections must be defined as: 



j^KLM ^ Fp{am,Yl-Fp{am,Y'' . R, = Zy + 0{g\. a^ 
I'-i/ciami, am2i 
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K2 


arueff 


Rq 


nKLM 

Kq 


Rp 


nKLM 
Up 


0.1150 


0.1150 


0.912 


0.999 


0.949 


1.026 


1.009 


0.1150 


0.1200 


0.822 


0.989 


0.942 


1.000 


0.991 


0.1150 


0.1250 


0.739 


0.981 


0.936 


0.978 


0.976 


0.1200 


0.1200 


0.731 


0.974 


0.934 


0.977 


0.974 


0.1200 


0.1250 


0.648 


0.962 


0.926 


0.958 


0.959 


0.1150 


0.1330 


0.618 


0.977 


0.929 


0.947 


0.950 


0.1250 


0.1250 


0.565 


0.946 


0.917 


0.937 


0.941 


0.1200 


0.1330 


0.528 


0.951 


0.916 


0.927 


0.932 


0.1150 


0.1425 


0.493 


0.996 


0.931 


0.919 


0.922 


0.1150 


0.1432 


0.484 


0.999 


0.932 


0.918 


0.920 


0.1250 


0.1330 


0.445 


0.927 


0.902 


0.910 


0.916 


0.1200 


0.1425 


0.402 


0.959 


0.913 


0.901 


0.903 


0.1200 


0.1432 


0.394 


0.962 


0.914 


0.899 


0.902 


0.1330 


0.1330 


0.324 


0.894 


0.881 


0.883 


0.889 


0.1250 


0.1425 


0.319 


0.924 


0.893 


0.882 


0.885 


0.1250 


0.1432 


0.310 


0.926 


0.894 


0.881 


0.883 


0.1330 


0.1425 


0.195 


0.872 


0.861 


0.855 


0.858 


0.1330 


0.1432 


0.190 


0.872 


0.860 


0.853 


0.855 


0.1425 


0.1425 


0.073 


0.823 


0.822 


0.823 


0.824 


0.1425 


0.1432 


0.065 


0.820 


0.819 


0.821 


0.821 


0.1432 


0.1432 


0.056 


0.817 


0.817 


0.818 


0.818 



Table 7: The ratios Rq and Rp, and their KLM-improved counterparts, for the Clover action 
at different values of the quark masses. The statistical errors are always less than (or equal 
to) 1 on the last digit. 



^ ryci{ami,am2) 
where the several KLM factors are given in eqs. (|100|), (|113|) and (120). Notice that, 



in the Clover case, the KLM corrections are always of 0(0^). 

The values of Rq, Rp and their KLM counterparts are given in table |^ for several 
values of the quark masses. The same ratios are also plotted in figures ^ and |^, 
with and without the KLM corrections, as a function of the effective mass am^ff = 
{ami + am2)/2. The main feature shown in figure ^ is that both Rq and Rp receive 
a large contribution from a term which is linear in the quark mass. This term affects 
the value of the two ratios by approximately 15% at am ~ 0.5 and by about 25% at 
am ~ 0.9. Moreover, such a large systematic error cannot be corrected by the KLM 
factors, the latter being at least quadratic in the quark masses. In fact, a qualitatively 
similar behaviour is found for the points in figure ^, where the KLM correction has 
been implemented. 

By looking at figure we also notice that the values of Rq, obtained with degenerate 
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Figure 3: The ratios Rq and Rp, for Clover fermions, as a function of the effective quark 
mass. We distinguish between the cases of degenerate and non- degenerate quarks in the 
current. 





ABC D 




0.80 0.19 -0.05 0.09 
0.80 0.16 -0.05 0.03 


Rp 

tjKLM 
tip 


0.81 0.12 -0.0008 -0.002 
0.80 0.14 -0.02 -0.008 



Table 8: Values of the coefficients A, B, C and D obtained from a fit of the ratios Rq, Rp 
and their KLM-improved counterparts. 



quark masses of the vector current, and those obtained by non-degenerate quarks, do 
not lie on the same curve. This suggests a significant dependence of this ratio on terms 
quadratic in the mass difference {ami — am2)^ . In this respect, figure ^ shows that 
the KLM correction improves the situation, since the difference between degenerate 
and non-degenerate points is largely reduced. However, the KLM prescription cannot 
correct large 0{g'^am) effects, which are evidently present in figure |^. It would also 
appear from this figure that such terms are independent of the ME's, since the difference 
between Rq and Rp settles only after masses of approximately am ~ 0.5. 

To make this discussion more quantitative, we fit the ratios Rq, Rp and their KLM 
counterparts, as in eq. ( [1251) . The resulting values of the coefficients A, B, C and D 
are given in table ^. The values of the coefficient A, in table ^, represent an estimate 
of Zv which is presumably free from any 0{am) corrections. We then quote: 

Zyff = 0.80(1) , from Rq , Rp (131) 
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Figure 4: Same as in figure 3 with KLM corrections. 



in remarkable agreement with its determination from light quark masses; cf. eq. ( |12g| ). 
Notice that the two results have been obtained from completely independent numerical 
simulations. 

The values of the coefficients B in table |8| provide an estimate of the size of the 
correction coming from a term which is linear in the quark mass. Since the action 
is tree- level improved, this term must be of 0{g'^a) or higher. We find B ~ 0.15, for 
both ratios Rq and Rp. As expected, these coefficients are practically unaffected by the 
KLM corrections. Remarkably, however, we find that their value is almost "universal", 
in the sense that it is almost independent on the external states of the ME of the vector 
current. 

On the other hand, the values of the coefficients C and D of the quadratic terms 
are not universal. For the ratio i?o, the large value of coefficient D is found to be 
significantly reduced by the KLM correction. 

We stress again that all these 0{a) effects are quite relevant in the region which, for 
current values of the lattice spacing, corresponds approximately to the charm quark 
mass. This in turn is crucial for the phenomenological studies of the heavy flavour 
physics on the lattice. 

In an attempt to reduce the large 0{g^a) effects observed in the data, we have also 
tried to modify the definition of the improved operators by varying the z parameter 
in eq. (84). Specifically, we have considered the values z = and z = 0.5 besides 
the standard choice z = 1 analyzed so far. We find that in all cases the size of the 
0{g'^a) contributions is approximately the same. Thus, no specific choice of z appears 
to improve the situation. 
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Ki 


A-2 


anieff 


Zv 




Zv (R^'^n 


0.1150 


0.1150 


0.912 


0.62(1) 


1.01(1) 


0.949 


0.1200 


0.1200 


0.731 


0.68(1) 


0.93(1) 


0.934 


0.1250 


0.1250 


0.565 


0.74(1) 


0.88(1) 


0.917 


0.1330 


0.1330 


0.324 


0.79(1) 


0.84(1) 


0.881 


0.1425 


0.1425 


0.073 


0.80(2) 


0.80(2) 


0.822 


0.1432 


0.1432 


0.056 


0.79(2) 


0.79(2) 


0.817 



Table 9: The RC Zy, for the Clover action, obtained from the axial WI with and without the 
KLM correction. For comparison, the values of the ratio Rq^^ are also shown in the last 
column. 



7.4 WI Determination of Zy and Za with the Clover Ac- 
tion and Heavy Quark Masses 

In this subsection we present the Clover results for the RC's Zy and Za obtained from 
a direct evaluation of the axial lattice chiral WI. This calculation has been performed 
by using the set of gauge configurations and quark propagators of subsection 

A determination of Zy can be obtained by considering the axial WI of eq. ([76|). 
In the Clover case, the right-hand side of this equation must be corrected for an 0{a) 
term as shown in eq. (^), so that the resulting WI has the form: 

2pi2 jdxj dy{P'\x)Af{y)vS\Q)) = "(^ " «p'') / dy{Af{y)A}?m (132) 

By computing p from eq. ( [6^ ) (with improved operators and symmetric derivative v'') 
and the two correlation functions appearing in the above WI, we can solve for Zy. 

The KLM improvement to eq. ( |132| ) involves both the standard KLM factors, which 
correct the two- and three-point correlation functions, and the KLM correction to the 
ratio yO, given by eq. ( |12lD . The results for Zy obtained from this WI are then presented 
in table with and without the KLM corrections. For a comparison, we also show in 
the table the values of Zy corresponding to the ratio Rq^^\ which have already been 
given in table |7[ 

Let us first consider in table ^ the values of Zy obtained from the WI (132) without 



the KLM corrections. As the bare quark mass increases from am ~ up to am ~ 0.9 
we find that Zy decreases by approximately 25% due to finite cutoff effects. This is to 
be contrasted to the increase by about 25% of the Rq estimate of Zy (sixth column of 
the table) in the same mass interval. The two determinations only agree at very small 
values of the quark mass. 

When we apply the KLM prescription we see from the fifth column of table ^ a 
dramatic change in the estimates of Zy. This effect is mainly due to the different 
normalization of the ratio p, and signals the presence of large 0{a'^m^) corrections. 
However, in spite of the KLM correction, the variation of Zy with the quark mass is 
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0.87(1) 


1.07(1) 
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0.1200 


0.731 


0.92(1) 


1.05(1) 


0.1250 


0.1250 


0.565 


0.97(1) 


1.05(1) 


0.1150 


0.1432 


0.484 


0.88(1) 


0.93(1) 


0.1200 


0.1432 


0.394 


0.91(1) 


0.95(1) 


0.1330 


0.1330 


0.324 


1.04(1) 


1.06(1) 


0.1250 


0.1432 


0.310 


0.95(1) 


0.97(1) 


0.1330 


0.1432 


0.190 


1.01(1) 


1.02(1) 


0.1425 


0.1425 


0.073 


1.11(2) 


1.11(2) 


0.1425 


0.1432 


0.065 


1.12(2) 


1.12(2) 


0.1432 


0.1432 


0.056 


1.12(2) 


1.12(2) 



Table 10: The RC Za for the Clover action, with and without KLM corrections, for different 
values of the quark masses. 



still of about 25%. Comparing the results of the last two columns of the table, we note 
a good agreement between the KLM corrected Zy estimates obtained with the WI 
and Ratio methods. This suggests that the KLM correction has adjusted the O(a^) 
terms in such a way that the resulting values of the RC's exhibit a more universal 
behaviour, approximately independent of the specific correlation functions used for 
their computation. 

Fitting Zv as in eq. ( |125D , we obtain an estimate of the RC's with presumably 
negligible 0{am) discretization effects. We quote: 

Z^ff = 0.80(2) , from WI (133) 

in perfect agreement with the previous determinations from the Ratio method (eqs. ( p.29D 
and (ll3ll)). 

The RC of the local axial current Za, can be computed from the WI of eq. (p8[), 
by choosing suitable spatial Lorentz indices of the currents: = p = 1,2,3. Also in 
this case we find it convenient to take the flavour -2 and -3 quarks degenerate in mass, 
so that the WI only involves the RC's Zj^ and Z'y. For the latter, we use the values 
obtained from the ratio Rq (at the corresponding quark masses), which is affected 
by a smaller statistical error, compared to the determination coming from the axial 
WI. Also for Za we have considered the results obtained with and without the KLM 
prescription, which involves in this case corrections to the correlation functions, to the 
ratio p and to Zy itself. The resulting values of Za are presented in table |l^ and are 
plotted in figure |5| as a function of the effective mass a^igyj = (ami + am2)/2. 

As in the case of the vector current, we find from table |l^ and figure ^ a large 
dependence of Za on the two quark masses. For values of the effective quark mass of the 
order of am ~ 0.5 such a dependence can affect the value of the RC by approximately 
20%. When we fit Za as in eq. (125), we find for the coefficient of the linear term 
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the value B ~ —0.13. This value is of opposite sign but similar in magnitude to the 
coefficient found for Zy- This indicates the presence of large 0{g'^am) corrections in 
both cases. From figure ^ we also notice that points corresponding to degenerate and 
non degenerate quark masses in the current do not lie on the same curve, indicating 
the presence of a large dependence on the difference of the two quark masses. From 
the fit we obtain D ~ —0.17 for both uncorrected and KLM corrected data. Therefore, 
we find that the KLM corrections do not improve the finite cutoff effects observed in 
the numerical estimates of Z^. 

From the zeroth order coefficient in the fit of Zj^ we also derive an estimate of the 
value of the RC with presumably negligible 0{am) effects. We quote: 

Z""/^ = 1.10(2) (134) 

in good agreement with the independent determination at light quark masses, presented 
in table ^. 

8 Further Improvements and Future Outlook 

The results presented in this paper lead to the following conclusions: 

• Applying various recipes of BPT gives RC's which are consistent to each other 
within 10% and which differ from those obtained from SPT by at most 20%. In 
most (but not all) cases, perturbative improvements based on BPT can lead to 
more reliable estimates of the RC's, in the sense that they are in better agreement 
with those obtained non-perturbatively. 
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At present values of the lattice cutoff, non-perturbative estimates of the RC's 
based on lattice WI's, suffer, in the Wilson case, from large discretization errors 
of about 20%. When tree-level Clover improvement is implemented, such sys- 
tematic effects decrease significantly for correlation functions computed at small 
quark mass. They remain large, however, as the quark mass increases {am ^ 0.2). 
We have checked that this situation persists for various choices of the field rota- 
tion, characteristic of tree- level 0{a) improvement (i.e. different choices of z in 
eq. 



• A very systematic analysis of the KLM corrections has been carried out, both 
with Wilson and Clover actions. In the Wilson case, the KLM factors sometimes 
succeed and sometimes fail to correct large 0{am) effects, depending on the 
correlation function examined. Moreover, even in the range of small quark masses, 
large residual 0{a) errors remain, as can be seen from the discrepancies in the Zy 
estimates from different correlation functions. In the Clover case, KLM cannot 
correct the observed large 0{gQam) systematic errors. 

The implication of these results is that, with present day lattice cutoffs of about 
2 — 4 GeV (/3 = 6.0 — 6.4), simulations of QCD phenomenology in the charm region 
(e.g. studies of or Z)-meson semileptonic decays) suffer from sizeable finite cutoff 
errors. This limits the predictability of lattice QCD for heavy flavour physics. In 
order to circumvent this problem, and given the impracticability of performing the 
simulations at signiflcantly smaller lattice spacings, one must seek to improve the 
lattice convergence to the continuum limit. In this work, we have examined two such 
efforts, namely Clover tree- level improvement and KLM matching. Here we comment 



briefly on the recently proposed extensions of these efforts, |32] and |39|. 

The aim of ref. is non-perturbative 0{a) improvement. To achieve this, the 
Wilson action must be modifled by adding the (dimension-5) Clover term of eq. ( [77[ ) 
with the coefficient csw determined non-perturbatively. In addition, a dimension-d 
renormalized operator must be cast in the form 

6 = Zoig\afi)[l + bo{g^)am][0 + aJ2cnig'')SOn] (135) 

n 

where the 50„'s are the dimension-(d -|- 1) operators which are allowed to mix with 
O. The constants bo and c„ can either be determined in PT or non-perturbatively. In 
PT they are known at tree-level and, in some cases, at 1-loop |5^. For example, for 
bilinear local quark operators, c„ = and bQ = 1 correspond to the tree-level Clover 
improvement discussed in this work, with field rotations determined at z = 1. For a full 
0{a) improvement, all the improving coefficients must be computed non-perturbatively. 



In ref. |32| the coefficients csw of the Clover term and ca_ of the axial current have been 
computed non-perturbatively, with the aid of lattice chiral WI's; for the coefficient cy 
of the vector current see ref. [^]. As shown in [^ ], chiral WI's cannot be used in the 
spirit of [^] for the determination of the coefficients bp (with the exception of by of 
the vector current computed in [pSl). Two methods for the computation of cq and bo, 
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which achieve 0{a) improvement at non-zero quark masses, can be found in refs 
and |35|. In this way, the approach to the continuum hmit is accelerated to O^a^). 



The proposal of ref. [^] can be summarized as follows. We saw that the KLM 
matching consists in the redefinition 

d = Zo{glafi)Fo{am)0 (136) 

The KLM factor Fq has been worked out at zero spatial momentum in the free theory 
and for any value of the quark mass. The authors of propose to extend this 
calculation beyond the free theory. They claim that it is possible to perform the lattice 
to continuum matching, again at zero spatial momentum but for any value of the quark 
mass, order by order in PT and to all orders in am. However, this matching has not 
yet been worked out beyond tree-level. 

We believe that the non-perturbative Clover improvement of on-shell ME's, pro- 



posed in 1 32], is the most promising way of controlling discretization errors at present 
day values of the UV cutoff. Our analysis of the tree-level Clover improvement at large 
quark masses strongly suggests that it is crucial to implement the extension of this pro- 



gram beyond the chiral limit |34, This will hopefully render lattice computations 



of heavy flavour physics extremely accurate. 
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Appendix: The Calculation of q* 

We give here the essential details and discuss some subtleties concerning the calculation 
of q* for the RC's Zy, Za and the ratio Zs/Zp. For definitiveness, we consider the 
MS renormalization prescription in the NDR scheme. The general 1-loop expression 
of eq. (p^), for the three quantities of interest, becomes 

Zv = 1 

Za = 1 

I- 



go(«)V . 

\2^P 



(47r) 
(47r)2 
(47r)2 



Cf 



r^LAT I r^LAT 



nLAT I riLAT 
^A + <-^S 



(137) 
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LAT 
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These expressions are easily derived from the discussion preceeding subsec.3.1. For 
convenience, we have rescaled the lattice constants C^^'^ by the colour factor Cf = 
4/3, which now appears explicitly. Note that none of the above quantities depends 
on the renormalization scale a/x. Consequently, the g*'s of these quantities are also 
independent of a^. 

For the Wilson action, the terms on the r.h.s of the above expressions have been 
worked out in ||2|, |^ : 



riLAT , r'LAT 



r^LAT I r'LAT 
^A + <-^S 

riLAT r'LAT 



1 1 



/3 + /s 



(138) 



-2h 



where the finite integrals /i, I2-, I3 and !■£ are reproduced here for completeness: 



2 r d\ 



/2 = 47rV2 



. (27r)4 



|(AiA4 - A5) + 2r^Ai + 4r^Af A4 ^ 4A4(A| - 2A1A3) 



Ij. = 47r 
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2 r 



AiAi 
4Ai(4- Ai) - A4 



A?Ai 



. (27r)4 

2 r 

-n (2^)4 



-4A 



3AiA2 



4AfA2 



^2 

5 , iA4 + 2r2A4 + 4r4Af 



+ 



(A3 



Ai 



,2 .2^ , A5 r^Ae Ai 
" + 4A?A2 " A2 " 2 



(139) 



with 



Ai =^sin2(g^/2) ; 
A5 = ^sin^iq/,) sm^{q^/2) 



A2 = ^sin2(g^)+4r2 
A4 = ^sin2(g^) 

; A6 = ^cos(g/, 



(140) 



In the present work, the Wilson parameter has been set to unity; r = 1. 

The recipe for finding q* is most straightforward for the ratio Zg/Zp. In order to 
implement eq. (|3l|), we must simply compute the 4-dimensional integral I2 and also 
the corresponding integral L2, obtained by multiplying the integrand of I2 by a \n{q^). 
L2 is also finite, and therefore calculable. Thus, for Zs/Zp we have 



(141) 
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We note in passing that the recipe of eq. (31) is not univocal, even for this simple case: 
the internal momentum q of the integral I2, being a dummy variable, can be chosen 
with a certain degree of arbitrariness (in particular, translations of q) which does not 
affect the final answer, since the integrand is made up of translationally invariant 
propagators. However, such changes in q will affect the logarithm in the integrand of 
L2, and thus will give us a different q* . This is only an academic problem, because q* 
is only meant to be a typical scale of the process under examination. Thus, estimates 
of q* which differ by the choice of internal momenta q are all equally acceptable. 

The calculation of q* for Zy and Za is more complicated. Besides L2, we also need 
Li and L3 (again obtained from Ii and I3, by modifying their integrands by ln(g-^)). 
Li and L3 are also finite, so this causes no serious problem. The real complication 
arises from the unity on the r.h.s of eq. ( |138| ). This constant is obtained as a difference 
of two integrals which diverge logarithmically with the external momentum ap^. We 
write them out explicitly. The first one arises from the divergent part of the projected 
amputated Green functions Tv{p) and Ta{p)- 



Iriap) 



167r^ 



^ d^q 



sm 



sm 



1 2 



sm 



2 (qu + apu 



V 2 

= - In(aV) --fE + i^oooi + 5/2 + 0{ap) 
whereas the second integral arises from the divergent part of the self-energy: 

d^q ^s.'^'^{qp + app)s\n{qp) 



(142) 



I^{ap) 



-IGvr^ 



-n (27r)^ 



4^ sinM^ 



sm 



qiy + apy 



In(aV) + 1E- i^oooi - 3/2 + 0{ap) 



(143) 



The integration (in the limit ap — > 0) has been performed with the aid of the tabulated 
integrals of |^8|. Thus, the factor +1 on the r.h.s. of the first two eqs. ( |138D is obtained 
from the integral I sum = It + It.'- 



I- 



SUM 



167r^ 



d'^q Ep [sill ilp) - sin (Qp + «Pp)] sin (gp) 



4E, 



sm 



1 2 r 



sm 



2 f qu + apu 



1 + 0{ap) 



(144) 



The above integral is finite. Its contribution to the RC's is given by its limit ap 0. 
We note that, in this limit the integrand of IsuM vanishes, except for q^ = 0. Thus, 
the integral is dominated by the region of small q. This enables us to substitute the 
integrand by its small q and ap limit and take it over the whole q range: 



ISUM 



1 



+00 



d\ 



ap • 



q^{q + apY 



(145) 
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We introduce the standard Feynman parameters, and after some trivial algebraic ma- 
nipulations obtain: 



IsUM 



— ^ / dx{l — x) 



ap-q 



[{q + xapY + {apYx{l — x)]' 



(146) 



We next shift the integration variable q ^ q + xap and pass over to polar coordinates 
in g-space. Integrating the polar angles leaves us with 



I SUM = '^{apf / dx{l 



+00 



dq^ 



(147) 



Finally, we perform a further change of variables 
out the integration in x to obtain 



[q"^ + {ap)'^x{l — x)]^ 

q^ /[a?p^x{l — x)], and carry 



IsUM 



dq' 



(g2 + 1)3 

The above is easily seen to be equal to 1, which is the final result of eq. (I43). 



(148) 



Now the application of eq. (31) requires the introduction of In(g^) in the integrand 
of I SUM- However, it is not obvious at which stage of the successive integrations from 
eq. ( |144| ) to eq. ( [148D this should be done. For example, if we were to introduce \n{q^) 
in eq. ( |147D , we would find a result proportional to In(a^p^). This does not make sense, 
since it yields a q* which depends on the momenta of the external states and, even 
worse, this ap dependence is singular. Thus we prefer to introduce the In(g^) in the 
integrand of eq. (148): 



LsUM 



+00 



dq 



2'f ln(g2) 



((?2 + 1)3 

The above integration can be easily performed; we find LsuM 
independent constant. Thus, q* for Zy is given by 



1 



ln(g 



* 2 \ 



1 1 



For Za we have 



ln(<i* 



4 2 



1 1 

-Li + -L2 
4 2 



(149) 

1, which is a p- 
(150) 



4 2 



/3 + /S 



(151) 



The Clover case presents no extra conceptual difficulties. Everything is like in the 
Wilson case, except that on the r.h.s. of eqs. ( |138| ) we have more finite integrals, which 
arise from the extra (Clover) term in the action and the Ip field rotations of subsection 
5.2 . These finite integrals are listed in refs. [0, 
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A final comment is in place here. A shortcut to all the subtleties related to LsuM 
would be to simply ignore the contribution of this integral to q* . In this case, q* is 
obtained from the last two equations, with the +l's omitted from both numerator and 
denominator. With this q* the most affected RC's are Zy and Za of the Clover action. 
Their values change by less than 2%, for all boosting recipes of the coupling. 
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